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A.  INTRODUCTION 

The  Final  Report  reports  on  the  activities  and  results  of  the  University  of  Illinois  for  the  US 
Army  Construction  Engineering  Research  Laboratory  Cond'act  Number  DACA8 8-94-0-0008, 
“Acoustical  Characterization  of  Soil.” 

Detection  and  classification  of  buried  cultural  artifacts  in  ground  soil  are  principal  gods  for 
the  production,  detection  and  processing  of  acoustic  signals.  Time  domain,  frequency  domain  and 
combined  time-frequency  domain  approaches  to  transmit  and^rocess  acoustic  signals  all  depend 
critically  on  the  acoustic  transduction  device  to  transmit  high-amplitude  acoustic  pressure  waves 
and  to  receive  low-amplitude  acoustic  pressure  waves  over  a  large  band  of  frequencies. 

The  transduction  device  is  the  means  to  an  end  and  represents  one  of  many  critical 
components  in  the  overall  system,  all  of  which  affect  the  system’s  full  utility.  Other  system 
components  include  the  transmit/rcceive  electronics,  the  broad-band  electrical  matching  between  the 
system’s  electronics  and  transduction  device,  the  low-power  and  high-power  switching  capability 
between  transmit  and  receive  when  the  same  transduction  device  is  used  for  transmit  and  receive, 
and  the  broad-band  mechanical  matching  between  the  transduction  device  and  the  acoustic 
propagation  medium. 

Once  a  cultural  or  archaeological  resource  site  is  identified,  it  must  then  be  assessed  in 
order  to  determine  its  significance  and  eligibility  for  National  Registry  of  Historic  Places 
(Executive  Order  1 1593).  A  Phase  II  eligibility  assessment  currently  costs  about  $10,000  to 
$30,000  per  site.  Given  that  there  are  approximately  120,000  archaeological  sites  in  the  Army 
alone,  the  cost  of  complete  Phase  n  assessments  is  prohibitive.  Therefore,  there  is  an  urgent  n^ 
to  significantly  reduce  the  cost  of  data  recovery  at  sites  with  an  unknown  probability  of  containing 
signSicant  cultural  or  archaeological  resources.  It  is  hypothesized  that  a  rapid  subsurface  acoustic 
imaging  technique  can  serve  this  need. 

The  investigators  wish  to  acknowledge  the  invaluable  assistance  of  Nadine  B.  Smith  and 
Richard  N.  Czerwinski  for  developing  the  data  acquisition  hardware  and  software,  Richard  N. 
Czerwinski  for  developing  the  data  analysis  software,  Catherine  Frazier  for  developing  much  of 
the  synthetic  aperture  image  approach.  Scott  Weisbrook  for  analyzing  the  soil  samples,  Dudley 
Swiney  for  andyzing  the  acoustic  data,  and  Scott  Weisbrook,  David  Tungate,  Dan  O’Brien,  Brett 
Boege  and  Kay  Raum  for  acoustic  data  acquisition. 

B.  MATERIALS  AND  METHODS 

1.  Acquisition  and  Evaluation  of  the  Soils 

Soils  samples  were  collected  from  Illinois  to  represent  a  wide  range  of  properties  expected 
to  influence  acoustic  response.  The  samples  were  chosen  to  include  low  and  high  contents  of 
organic  matter,  sand,  silt  and  clay,  and  a  range  of  clay  mineralogy.  Table  1  gives  the  locations, 
depths,  horizons,  and  classifications  of  the  soils. 

Pedologists  recognize  layers  in  soils  that  are  called  horizons.  These  can  be  of  various 
depths  and  thicknesses,  but  the  properties  of  the  soil  material  within  individual  horizons  are 
different  from  the  other  horizons  in  a  soil  but  internally  similar.  Horizons  are  labeled  “A”,  “B”, 
and  “C”  down  from  the  surface.  The  horizon  labels  are  given  subordinate  distinctions  to  indicate  a 
particular  subtype  of  horizon.  The  term  “Ap”  refers  to  the  surface  horizon  of  a  soil  that  has  been 
plowed  (the  “p”  designation)  or  otherwise  disturbed  sometime  in  the  past  Most  surfaces  of  soils 
have  been  plowed  or  disturb^  by  forestry  or  agriculture  to  the  point  that  soil  scientists  recognize 
the  surface  horizon  (the  “A”)  as  an  Ap  and  not  simply  as  an  “A”  horizon  that  would  indicate  that  it 
was  “natural”  or  uncultivated  or  undisturbed  by  human  activities. 
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Table  1.  Locations  and  classifications  of  soils  used  for  acoustic  characterization. 


Sample 

Code 

Soil 

Series 

Soil 

Classification 

Soil 

Horizon 

Collection 
Depth  (cm) 

HSII 

1 

CAB 

Catlin 

fine-silty,  mixed, 
mesic. 

Tvpic  Argiudoll 

Bt 

73-85 

Champaign 

40  05  14 

88  13  55 

2 

PLA 

Plain- 

field 

mixed,  mesic. 
Typic  Udipsamment 

Ap 

0-20 

Mason 

40  09  28 

90  05  14 

3 

SAC 

Sable 

Cg 

300-320 

hfcnaid 

40  00  00 

88  42  30 

■ 

DRA 

Dnim- 

mer 

fine-silty,  mixed, 

mesic. 

Typic  Haplaquoll 

Ap 

0-20 

Champaign 

40  05  01 

88  13  55 

5 

ADA 

Adrian 

sandy,  mixed,  mesic, 
euic. 

Terric  Medisaprist 

■ 

0-20 

Cass 

40  0100 

90  24  00 

6 

MEA 

Med¬ 

way 

fine-loamy,  mbced, 
mesic.  Fluyaquentic 
Hapludoll 

Ap 

0-20 

Cass 

40  00  00 

90  28  00 

Prior  to  analysis,  the  soil  samples  were  air  dried  and  sieved  to  exclude  the  >2  mm  fraction. 
The  influence  of  coarse  fragment  content,  structure,  and  consistency  was  not  addressed  by  this 
phase  of  the  work.  The  soil  samples  were  characterized  as  to  particle-size  distribution,  organic 
matter  content,  clay  mineralogy,  coefficient  of  linear  extensibility,  plastic  and  liquid  limits, 
plasticity  index,  and  electrical  conductivity.  Table  2  gives  the  soil  physical  characterization  data. 


Table  2.  Physical  properties  of  acoustically  characterized  soils. 


Soil 

Code 

Sand 

% 

Silt 

% 

Qay 

% 

Organic 

Matter 

% 

Soil 

Texture  Class 

Liquid 

limit 

Plastic 

limit 

Plasticity 

index 

Conductiyity 

dS/m 

COLE 

% 

CAB 

11 

El 

36 

K9i 

43 

24 

19 

143 

14 

PLA 

97 

1 

2 

0.4 

Sand 

NP 

NP 

NP 

120 

NP 

SAC 

2 

82 

16 

0.1 

Silt  Loam 

29 

24 

5 

379 

4 

DRA 

12 

50 

38 

9.8 

Silty  Clay  Loam 

43 

28 

15 

363 

13 

ADA 

72 

18 

10 

11.7 

Fine  Sandy  Loam 

NP 

NP 

NP 

688 

NP 

MEA 

38 

38 

24 

2.3 

Loam 

28 

19 

9 

258 

7 

Soil  characterization  was  by  standard  methodology.  Particle-size  analysis  was  determined 
by  the  hydrometer  method  (Gee  and  Bauder,  1986).  The  silt  and  clay  were  determined  by 
hydrometer  and  the  sand  fraction  was  separated  by  sieving.  Electric^  conductivity  of  the  soils  was 
determined  on  a  1:1  extract  (Dahnke,  1988). 

Organic  carbon  content  was  determined  by  wet  combustion  with  sodium  dichromate 
(Nelson  and  Sommers,  1986).  Organic  matter  was  calculated  from  the  organic  carbon  content  by 
multiplying  by  a  conversion  factor  of  1.724.  Liquid  and  plastic  limits  and  plasticity  index  was 
determined  with  a  glass  plate  and  an  electric  liquid  limit  machine  (AASHO,  1969).  Coefficient  of 
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linear  extensibility  (COLE)  was  determined  by  air  drying  wet  samples  created  with  an  extrusion 
device  (Schafer  and  Singer,  1976). 

Mineralogy  is  being  determined  on  the  sUt  and  clay-sized  material  by  X-ray  diffi^ction 
(Whittig  and  Allardice,  1986).  The  size  fractions  were  separated  by  centrifugation  and  analyzed  on 
Ae  Department  of  Natural  Resources  and  Environmental  Sciences’  X-ray  diffractometer. 

Water  content  of  the  soils  was  determined  by  weight  loss  on  drying,  then  conversion  to 
volumetric  water  content  basis  by  multiplying  by  the  bulk  density  (Ganlner,  1986). 

The  original  experimental  design  detailed  six  soils,  three  compaction  levels,  and  four 
moisture  levels  for  a  total  of  72  experimental  units.  However,  the  Plainfield  soil  was  evaluated 
only  loosely  compacted  at  the  low  moisture  level,  and  experimental  difficulties  prevented  three 
levels  of  compaction.  Therefore,  only  the  two  extreme  compaction  levels  were  evaluated,  viz., 
loose  and  dense.  Also,  for  the  100%  moisture  condition,  it  was  necessary  to  measure  only  one 
compaction  level  because  of  the  incompressibility  of  the  water  that  saturated  the  soil  pores. 
Therefore,  the  final  experimental  design  included  five  soils  at  two  compaction  levels  and  three 
moisture  levels,  plus  dry,  loose  Plainfield,  plus  Plainfield  at  two  moistures  and  two  compactions, 
plus  all  six  saturated  soils.  This  gave  a  total  of  41  experimental  units.  The  soils  were  acoustically 
evaluated  at  a  minimum  of  three  thicknesses  ranging  from  5-27  cm.  Soil  thickness  was  varied  both 
by  addition  of  new  soil  to  the  container  and  by  removal  of  material.  The  total  number  of  acoustic 
evaluations  was  231  (Appendix  Al). 

The  soil  was  placed  in  a  calibrated  plastic  tub  (approx.  30  cm  high  by  30  cm  diam.)  and 
weighed  prior  to  acoustic  evaluation.  A  small  sample  was  removed  to  determine  gravimetric 
moisture  content  at  the  time  the  acoustic  evaluation  was  made.  Soil  particle  density  was  estimated 
by  adjusting  the  density  of  quartz,  2.65  g/cc,  by  the  soil  organic  matter  content  at  a  density  of  1.3 
g/cc.  From  these  data  the  soil  sample  bulk  density,  moisture  content,  and  thickness  were 
determined. 

Bulk  density  of  the  soils  was  determined  by  weight  per  unit  volume  (Blake  and  Hartge, 
1986).  Porosity  of  the  soils  was  determined  by  calculating  the  ratio  of  sample  weight  to  volume 
(Danielson  and  Sutherland,  1986). 

Initial  water  content  of  the  soils  for  acoustic  evaluation  was  air-dry.  The  two  next  higher 
water  contents  were  nominally  increased  to  10  and  25%  water  by  volume.  This  was  achiev^  by 
carefully  adding  the  appropriate  quantity  of  water  and  thoroughly  mixing  and  screening  the  moist 
soil  through  a  4  mm  sieve.  The  saturate  soil  was  made  by  adding  screened  moist  soil  to  standing 
water.  All  moist  soils  were  allowed  to  equilibrate  for  at  least  24  hours  before  acoustic  evaluation. 

The  loose  compaction  of  the  soils  was  achieved  by  pouring  the  sieved  soil  into  the 
calibrated  tub  and  striking  off  a  level,  smooth  surface.  The  soils  when  air-dry  were  quite  powdery 
and  dusty.  Compaction  was  achieved  by  adding  a  succession  of  2  cm  lifts  to  the  bottom  of  the  tub 
and  applying  a  load  by  means  of  vigorous  stomping  onto  a  plywood  disk  on  the  soil  surface.  This 
was  repeated  until  the  desired  soil  thickness  was  achieved. 

The  acoustic  data  was  recorded  in  files  that  were  named  to  identify  the  experimental  soil 
conditions.  The  first  three  letters  in  the  file  name  were  the  soil  code  (see  Appendix  Al).  The 
next  two  numbers  were  the  sample  thickness  in  cm.  The  next  digit  was  the  moisture  code,  1, 2,  3, 
or  5  for  air-dry,  10%,  25%,  and  saturated,  respectively.  The  next  digit  was  the  compaction  code, 

1  for  loose,  4  for  compact.  The  last  number  in  the  file  name  represented  the  replication  of  those 
soil  conditions.  An  additional  code  is  added  after  the  file  name  to  indicate  sequence  number.  Table 
Al  lists  the  sample  files  and  the  soil  conditions  at  the  time  the  acoustic  data  was  collected. 
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A  block  diagram  of  the  soil  characterization  system  is  shown  in  Figure  1.  The  system 
consists  of  a  host  computer  (ZEOS  133  MHz  Pentium,  Minneapolis,  MN)  which  is  used  to  control 
all  data  acquisition  procedures.  This  computer  has  a  16  MB  of  core  memory  (expandable  to  32 
megabytes),  700  MB  hard  drive^torage,  Mitsumi  4X  CD-ROM  Drive,  SuperVGA  graphics,  and  is 
connected  to  the  Beckman  Institute  internet  network.  Communications  witii  the  internet  network  is 
through  NCSA  Telnet  2.6.1.  A  Hewett  Packard  HP81 16A  (with  option  001)  programmable  signal 
generator  receives  instructions  via  the  IEEE-488  communication  bus  which  prcxiuces  a  low-level 
signal.  The  low-level  signal  from  the  HP8 1 16A  is  amplified  by  a  3000  W  power  amplifier 
(Industrial  Test  Equipment  Powertron  3000S  amplifier.  Port  Washington,  NY)  and  impedance 
matched  to  the  N^  F56  Serial  58  acoustic  source  (Transducer  Branch,  U.S.  Naval  Research 
Laboratory’s  Underwater  Sound  Reference  Detachment,  Orlando,  FL)  by  a  transformer  (Industrial 
Test  Equipment  ET-900V  transformer). 

The  acoustic  signal  is  propagated  through  the  soil  sample  contained  in  a  tub  which  is 
coupled  to  face  of  the  F56  acoustic  source  via  a  water  interface.  The  transmitted  acoustic  signal  is 
received  by  the  NRL  F42C  Serial  28  hydrophone  (Transducer  Branch,  U.S.  Naval  Research 
Laboratory’s  Underwater  Sound  Reference  Detachment,  Orlando,  FL)  which  is  acoustically 
coupled  to  the  top  of  the  soil  sample  with  DOW  7 10  pheynlated  silicon  oil.  The  NRL  F42C 
hydrophone  is  stable  and  calibrated  with  a  bandwidth  between  100  Hz  and  200  kHz.  The 
c^ibration  is  traceable  to  the  NRL’s  Underwater  Sound  Reference  Detachment,  a  national 
standards  laboratory. 

The  backscattered  acoustic  signal  is  acquired  by  two  SE1025-H  acoustic  emission 
hydrophones  (Dunegan  Engineering,  San  Juan  Capistrano,  CA). 

The  four  low  level  signals  from  the  HP  81 16A  signal  generator  (Hewlett-Packard),  the 
NRL  F42C  hydrophone  and  two  SE1025-H  acoustics  emission  hydrophones  are  amplified  with 
gains  of  10, 100, 500  and  500,  respectively.  The  design  of  the  high  input  impedance  amplifiers  is 
shown  in  Figure  2.  The  outputs  from  the  amplifiers  are  digitized  by  a  DT2812A  A/D  board  (Data 
Translation,  Marlboro,  MA). 

Data  Acquisition  Software 

A  set  of  programs  for  the  data  acquisition  system  instructs  the  signal  generator  to  scan  a 
range  of  amplimdes  and  frequencies  and  to  receive,  digitize  and  store  in  individual  ASCII  files  the 
four  low  level  signals. 

Data  acquisition  and  command  instructions  to  the  IEEE-488  compatible  devices  from  the 
ZEOS  Panthera  occurs  thorough  the  General  Purpose  Interface  Bus  (GPIB)  (National  Instruments- 
GPIB/PC,  1989).  The  interface  bus  is  comprised  of  a  interface  board  (National  Instruments  PCII) 
installed  in  the  motherboard  of  the  ZEOS  computer  and  physically  connected  through  a  GPIB  cable 
to  the  electronic  devices  which  allows  the  electronic  devices  to  communicate.  A  program 
GPIB.COM  (supplied  by  National  Instruments)  is  a  software  loadable  device  driver  which  is 
loaded  at  the  system  start-up  by  DOS. 

All  of  the  communication  programs  and  the  autoexec.bat  file  (Appendix  A2)  are 
specifically  written  for  use  of  the  ^crosoft*  C6.0.  Within  the  autoexec.bat  file  are  instructions 
to  the  operating  system  as  to  the  location  of  libraries  and  include  files.  The  config.sys  file 
(Appendix  A3)  is  also  specially  written  for  the  computer  environment. 
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10  kOhms 


Vcc=+14V 


LM351 


To  the  Data  Translation 
DT2812A  Board 


100  kOhms 


220kOhms 


Vcc=-14V 


Amplified 

Device 

HP8116A 

NRL42C 

SE1025H-1 

SE1025H-2 


Amplifier  No# 
0 


Rf  Value 

loom 

1.0  MQ 
4.7  MQ 
4.7  Ma 


DT2818A 
Channel  No# 
0 


Figure  2 

A  schematic  of  the  high  input  impedance  amplifier  circuit  for  the  HP8 1 1 6A  signal  generator 
source  and  the  three  hydrophones  G'JRL  42C,  SE1025H-1  and  SE1025H-2).  The  table 
lists  the  desired  gain  for  the  individual  device  and  the  require  feedback  resistance  value  for 
Rf.  The  table  also  indicates  the  input  channel  for  the  for  the  Data  Translation  DT2812A 
A/D  board. 
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The  soil4.mak  file  (Appendix  A4)  is  the  make  file  for  the  soil4.exe  program.  The 
puipose  of  the  make  file  is  to  compile  all  the  source  code  for  the  program,  define  the  type  of 
memory  model  to  use,  create  the  stack  of  32()(X)  bytes  and  then  link  all  of  the  object  files. 

For  use  of  the  software  which  runs  the  soil  characterization  system  are  3  non  Microsoft* 
C6.0  include  files.  The  DTSTJTHS.H  program,  written  by  Data  Translation  (Appendix  A5), 
contains  a  list  of  define  (#define  in  C  language)  for  the  Data  Translation  software  tools  called  fiom 
a  Microsoft*  C6.0  program.  The  DTST_XMALH  (Appendix  A6)  program  is  the  definition  file 
for  the  extended  memory  management  functions  using  HIMEM.SYS  from  Microsoft*  C6.0. 
Finally,  the  DTSTJERR.H  (Appendix  A7)  program  contains  the  DTI  Software  Tools  error 
codes  which  provides  4096  different  error  codes. 

Data  translation  supplies  two  libraries  which  contain  error  information  and  codes  but  which 
are  too  extensive  to  list  in  the  appendix.  The  C  program  which  calls  these  libraries  is  called 
MSC_SUB.C  (Appendix  A8).  This  program  contains  the  source  which  basically  creates, 
calls,  clears,  utilizes  and  deletes  the  buffers  which  store  the  acquired  data  from  the  channels. 
Detailed  information  about  using  the  Data  Translation  libraries,  include  and  source  code  can  be 
found  elsewhere  (Data  Translation  SPO  126  Software  Tool  Kit,  1992). 

The  soil4.c  (Appendix  A9)  program  is  the  source  code  for  the  executable  main  program 
which  controls  the  soil  ciata  acquisition. 

The  procedure  to  transfer  raw  data  from  soil  acquisition  system  (bmode.brl.uiuc.edu)  to 
workstation  (ecstasy)  is  listed  in  Appendix  AlO. 

Data  Processing  Software 

The  processing  of  the  data  is  done  in  Matlab  (The  MathWorks,  Natick,  MA)  on  a 
SUNSparc2  workstation.  The  analysis  code  is  modular  in  design  to  facilitate  experimental 
analysis.  Processing  is  performed  by  a  succession  of  two  primary  routines,  procsoil.m  and 
process. m. 

The  procsoil.m  (Appendix  All)  program  reads  in  (v/a  loadsoil.m.  Appendix 
A 12)  the  raw  data  soil  measurements  from  a  single  thickness,  and  produces  a  list  of  estimates  of 
the  time  delay  between  pulse  transmission  and  reception,  and  the  amplitude  of  the  received  pulse. 
The  estimates  are  tabulated  in  1  kHz  increments  of  frequency,  and  each  value  is  obtained  from  four 
different  frequency  measurements,  i.e.,  the  1  kHz  values  are  obtained  from  the  scans  at  1000  Hz, 
1250  Hz,  1500  Hz  and  1750  Hz,  from  1  kHz  to  10  kHz. 

Time  of  flight  (TOF)  is  estimated  by  computing  the  correlation  function  between  the 
transmitted  and  received  pulses  for  each  of  the  frequencies  being  combined,  and  forming  a 
histogram  of  the  peak  locations  which  occur.  Two  estimates  represent  high  and  low  values 
intended  to  define  a  range  of  possible  TOF  values,  a  third  estimate  is  the  median  of  the  first  two, 
and  the  fourth  estimate  is  the  most  frequently  occurring  histogram  peak  location.  In  practice,  this 
should  be  the  most  accurate  estimate. 

The  amplitude  values  returned  by  procsoiLm  are  correlation  function  values 
corresponding  to  the  TOF  estimates.  These  values  are  directly  proportional  to  the  desired  quantity, 
with  a  scaling  factor  equal  to  the  square  of  the  transmitted  pulse  amplitude;  as  long  as  data  sets  of 
different  transmitting  signal  amplitude  are  not  mixed,  there  is  no  ne^  to  normalize  measurements 
by  this  quantity. 

The  peak.m  (Appendix  A13)  program  is  a  function  which  does  the  work  on  the  signals 
within  procsoil.m  (except  for  the  portions  dealing  with  user  interface).  It  allows  for  greater  ease 
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of  software  development  by  containing  the  commands  frequently  invoked  in  procsoil.m  into  a 
single  file. 

The  process.m  (Appendix  A14)  program  computes  frequency  dependent  speed  of 
sound  estimates  by  performing  linear  regressions  with  respect  to  thickness.  Speed  of  sound  is 
calculated  as  the  slope  of  line  which  best  fits  (in  a  least  squares  sense)  the  plot  of  TOP  vs. 
thickness.  Attenuation  is  the  slope  of  the  line  which  best  fits  the  plot  of  20*log(received 
amplitude)  vs.  thickness.  The  20  logO  operation  is  performed  so  that  attenuation  can  be  reported 
in  &/cm-kHz.  The  process.m  program  also  reports  the  coefficient  of  the  regression  to  indicate 
how  well  the  points  fit  the  regression  line. 

The  process.m  program  also  reports  the  frequency  dependent  results  from  the  log 
decrement  analysis.  The  log_dec.m  (Appendix  A15)  program  is  a  function  used  to  analyze  the 
data  received  by  the  acoustic  emissions  transducers,  positioned  at  the  bottom  of  the  mb  of  soil.  It 
computes  the  center  frequency  and  3  dB  bandwidths  of  the  exponentially  decaying  signal  which 
appears  after  the  excitation  signal  is  turned  off.  By  the  Biot  theory  of  wave  propagation  in  porous 
m^ia,  the  3  dB  bandwidth  is  related  to  the  decay  rate  of  the  decaying  component;  this  decay  rate  is 
calculated  direcdy  and  ft'om  the  3  dB  bandwidth. 

The  coef.m  (Appendix  A16)  program  is  a  function  which  evaluates  the  correlation 
function  between  two  data  sequences  at  a  certain  point.  The  maxes.m  (Appendix  A17) 
program  is  a  function  which  return  the  indices  of  local  maxima  in  a  data  sequence.  It  is  used  here 
to  determine  the  positions  of  correlation  peaks.  The  regression.m  (Appendix  A18)  program  is 
a  function  used  to  compute  the  best  fit  line,  in  a  least  squares  sense,  to  a  set  of  points. 

3.  Acoustic  Measurement  Technique 

Samples  were  hand-packed  in  the  soil  containment  vessel  for  acoustic  evaluation.  This 
hand-packed  sample  was  at  a  relatively  low  density  (loose  compaction).  After  this  sample  was 
evaluated  acoustically,  the  sample  density  (compaction)  was  increased  by  placing  a  piston-type 
aixangement  on  top  of  the  soil  in  the  soil  containment  vessel  and  applying  an  impulsive-type  force 
to  the  piston. 


C.  SOIL  RESULTS  AND  DISCUSSION 

Tables  3  and  4  provide  the  mean  values  of  attenuation  coefficient  (dB/cm-kHz)  and 
propagation  speed  (m/s)  for  the  six  soil  types  (ADA,  CAB,  DRA,  MEA,  PLA,  SAC)  as  a  function 
of  four  soil  moistures  and  two  soil  compactions.  The  first  digit  of  the  code  in  the  left  hand  column 
designated  soil  moisture  (1, 2,  3, 5)  where  1  designated  dry  soil  and  5  designates  fully 


Table  3.  Mean  value  of  the  attenuation  coefficient  (dB/cm-kHz) 


ADA 

CAB 

DRA 

MEA 

PLA 

SAC 

11 

0.31 

0.14 

0.21 

0.12 

0.23 

0.39 

14 

0.20 

0.59 

0.53 

0.63 

0.36 

21 

0.68 

0.36 

0.38 

0.39 

0.35 

0.46 

24 

0.58 

0.91 

0.56 

0.62 

0.50 

0.69 

31 

0.39 

0.45 

0.51 

0.28 

0.81 

0.50 

34 

0.71 

0.90 

0.71 

0.90 

0.16 

0.96 

51 

0.49 

0.58 

0.58 

0.31 

0.75 

0.38 
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Table  4.  Mean  value  of  the  propagation  speed  (m/s) 


ADA 

CAB 

DRA 

MEA 

PLA 

SAC 

11 

153 

190 

177 

154 

260 

140 

14 

121 

150 

159 

126 

— 

139 

21 

89 

114 

118 

151 

138 

117 

24 

147 

102 

136 

93 

103 

122 

31 

87 

150 

107 

161 

122 

121 

34 

86 

105 

245 

105 

253 

176 

51 

102 

102 

118 

139 

189 

207 

saturated  soil.  The  second  digit  of  the  code  in  the  left  hand  column  designates  soil  compaction 
(1,4)  where  1  designates  loose  compaction  and  4  designates  dense  compaction.  Appendix  A19 
provides  a  more  complete  listing  of  the  acoustic  properties  as  a  function  of  the  detailed  soil 
characteristics. 


Figures  3-9  graphically  represent  the  attenuation  coefficient  (dB/cm-kHz)  and 
propagation  speed  (m/s)  for  each  of  the  seven  soil  conditions  (11, 14,  21, 24,  31,  34,  51).  The 
last  digit  in  the  code  in  the  figures  is  a  sequence  number  which  was  used  for  repeat  experimental 
units. 


The  experimental  protocol  called  for  measuring  the  acoustical  properties  in  relatively 
homogeneous  soil  samples.  It  is  thus  suggest  that  the  experimental  approach  for  homogeneous 
soil  preparation  is  the  principal  factor  which  contributed  to  the  relatively  narrow  range  of  measured 
acoustic  propagation  propeny  values.  Additionally,  because  the  soil  samples  were  relatively 
homogeneous,  this  precluded  an  evaluation  of  the  acoustic  scattering  properties  because  scattering 
results  from  acoustic  heterogeneities. 

The  attenuation  coefficient  over  the  1-10  kHz  frequency  range  ranged  between  a  low  of 
about  0.1  dB/cm-kHz  which  was  most  prevalent  for  the  loose  dry  (11)  samples  and  a  high  close  to 
1  dB/cm-kHz  which  was  most  prevalent  for  the  more  moist  compact  (34)  samples. 

The  attenuation  coefficient  has  a  direct  influence  over  the  imaging  depth  for  a  given 
dynamic  range.  Figure  10  demonstrated  the  influence  of  the  roundtrip  propagation  loss  as  a 
function  of  attenuation  coefficient  (top  panel)  and  frequency  (bottom  panel).  For  a  dynamic  range 
of  140  dB  (typical  of  diagnostic  ultrasound  imaging  systems  but  unknown  yet  for  a  sub-surface 
acoustic  imaging  system),  the  top  panel  in  Figure  10  suggests  an  imaging  depth  of  4.6  feet  is 
achievable  for  an  attenuation  coefficient  of  0.5  dB/cm-kHz  at  a  frequency  of  1  kHz  and  an  imaging 
depth  of  6  feet  is  easily  achievable  for  lesser  attenuation  coefficient  values.  For  the  same  dynamic 
range,  the  bottom  panel  in  Figure  10  suggests  an  imaging  depth  of  2.3  feet  is  achievable  for  a 
frequency  of  5  kHz  at  an  attenuation  coefficient  of  0.2  dB/cm-kHz,  increases  to  3.8  feet  for  a 
frequency  of  3  kHz  and  to  more  than  6  feet  for  a  frequency  of  1  kHz. 

It  is  also  necessary  to  consider  the  acoustic  reflection  coefficient  from  the  object  to  be 
imaged  in  soil  when  considering  the  overall  roundtrip  propagation  loss,  but,  in  general,  it  appears 
that  the  acoustic  reflection  coefficient  will  be  around  0  dB.  TTiis  is  based  on  the  experimental 
results  obtained  herein.  The  characteristic  acoustic  impedance  of  soil  is  in  the  range  of  (1-3)  x  10* 
Pa-s/m  (density  «10(X)  kg/m*;  propagation  speed  » 100-300  m/s).  A  plastic  object  might  have  a 
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characteristic  acoustic  impedance  in  the  range  of  3.2  x  10*  Pa-s/m  (for  Lucite;  density  =  1200 
kg/m^;  propagation  speed  =  2650  m/s)  whereas  a  metallic  object  would  have  a  greater  value  of 
characteristic  acoustic  impedance.  The  pressure  reflection  coefficient,  at  norm^  incidence,  is 

7  —  7 

—  otiject  ■^soil 
^object  ^soil 

where,  for  =  3.2  x  10*  Pa-s/m  and  =  3  x  10*  Pa-s/m,  R  =  0.826  (or  -1.6  dB). 

The  imaging  depth  is,  thus,  inversely  proportional  to  frequency  which  results  in  the 
important  engineering  tradeoff  between  depth  of  penetration  into  soil  and  image  resolution.  As  a 
crude  rule  of  thumb  (dealt  with  more  precisely  in  Section  D),  the  image  resolution  can  be 
approximated  to  that  of  the  acoustic  wavelength  ( A  =  c/f )  where  c  is  the  propagation  speed  and  f 
is  the  acoustic  frequency.  Thus,  as  frequency  increases,  the  wavelength  decreases  (resolution 
improves)  and  depth  of  penetration  decreases. 

A  detailed  analysis  to  the  mechanisms  responsible  for  the  interaction  of  the  propagated 
acoustic  wave  with  soil  is  beyond  the  scope  of  the  project  However,  the  initial  hypothesis 
suggested  that  the  acoustic  propagation  properties  of  in  soil  might  be  a  function  of  soil  moisture 
and  compaction.  A  first  order  examination  of  the  acoustic  propagation  properties  as  a  function  of 
moisture  and  compaction  (see  Figures  11-16)  do  not  reveal  any  obvious  trends. 

D,  SUB-SURFACE  IMAGING 

1.  Analytical  component 

The  subsurface  image  formation  problem  has  multiple  strategies.  Whether  performing 
conventional  B-scan  or  using  synthetic  apenure  techniques  to  create  an  image,  the  data  are 
collected  by  scanning  the  transducer  or  the  transducer’s  beam  along  a  line.  For  B-scan  imaging, 
the  positions  where  the  RF  data  are  collected  must  be  far  enough  apart  that  the  signals  can  be 
considered  independent  of  each  other  while  still  satisfying  the  spatial  sampling  criterion.  Each 
returned  (backscattered)  acoustic  echo  is  detected  and  the  image  is  created  by  stacking  the  detected 
scans  side  by  side  so  that  one  of  the  scan  dimensions  represents  lateral  position  of  the  scan,  and  the 
other  represents  the  axial  depth  into  the  medium,  and  the  pixel  brighmess  represents  the  amplitude 
of  the  detected  acoustic  echo.  In  this  case,  the  ideal  beam  would  be  very  narrow  to  yield  good 
lateral  resolution,  and  the  ideal  transmitted  pulse  would  have  narrow  time  duration  to  yield  good 
axial  resolution.  Within  certain  limits  of  approximation,  lateral  resolution  is  that  of  the  beamwidth 
and  described  mathematically  at  the  focus  as 


Lateral  Resolution  = 


A  =  f/#A 


(2) 


where  FL  is  the  focal  length,  D  is  the  source  diameter,  A  is  the  acoustic  wavelength  and  fr#  is  the 
f-number  (ratio  of  focal  length  to  source  diameter).  In  practical  situations,  the  be^  is  narrow  only 
over  a  small  depth  near  its  focus.  This  means  that  only  part  of  the  image  is  actually  in  focus.  To 
compensate  for  this,  the  received  echoes  can  be  dynamically  focused  on  receive,  thus  creating 
multiple  focal  regions,  each  described  by  a  different  lateral  resolution.  Figure  17  schematically 
shows  three  focal  regions  for  the  same  source.  The  final  image  results  from  a  composite  of 
multiple  focal  regions,  thus  optimizing  (minimizing)  the  lateral  resolution  along  a  single  B-scan 
image  line. 
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Also  within  certain  limits  of  approximation,  axial  resolution  is  that  of  the  space  occupied  by  an 
acoustic  pulse  and  described  mathematically  as 


Axial  Resolution  =  QA  = 


c 


(3) 


where  the  system  Q  is  the  ratio  of  center  frequency  f  to  system  bandwidth  Af  and  c  is  the  acoustic 
propagation  speed. 

In  the  synthetic  aperture  case,  the  RF  data  are  collected  using  a  transducer  with  a  wide 
beam  pattern  and  a  correspondingly  small  aperture.  Again,  the  ideal  pulse  has  a  narrow  time 
duration  (or  large  bandwidth).  RF  data  are  collected  at  lateral  positions  close  enough  together  that 
the  area  of  interest  is  illuminated  by  many  transmit  beams.  In  order  to  synthesize  the  wide  aperture 
of  a  focused  transducer,  several  neighboring  beams  are  coherently  summed.  Obtaining  a  focused 
beam  at  a  given  position  is  a  matter  of  making  sure  the  beams  add  constructively  at  that  point  or,  in 
other  words,  have  the  same  phase  at  that  point  This  can  be  accomplished  by  delaying  the  beams 
near  the  center  of  the  sampled  aperture.  It  is  proposed  to  adapt  signal  processing  teclmiques  from 
synthetic  aperture  radar  (SAR). 

SAR  is  a  microwave  imaging  technique  that  produces  imagery  ha\ing  better  lateral 
resolution  than  that  dictated  by  the  source’s  beamwidth.  In  the  most  widely  known  form  of  SAR, 
called  strip-mapping,  the  source’s  beam  sweeps  out  a  strip  on  the  ground,  a  sequence  of  closely 
spaced  pulses  is  transmitted,  and  the  returned  waveforms  are  recorfed  (Munson  and  Visenten, 
1989).  The  source  moves  only  a  small  fraction  of  its  beamwidth  between  pulse  transmissions  so 
that  each  point  in  the  scene  is  illuminated  by  many  pulses.  The  returned  RF  signals  are  correlated 
in  the  range  direction  with  a  replica  of  the  transmitted  signal,  and  in  the  cross-range  direction  with  a 
linear  FM  waveform  to  produce  a  high-resolution  radar  image.  Within  certain  limits  of 
approximation,  utilizing  this  form  of  processing  with  a  source  diameter  D  yields  lateral  resolution 
of  D/2.  Thus,  a  smaller  source  actually  provides  BETTER  lateral  resolution!  It  is  our  hypothesis 
that  a  similar  form  of  processing  can  be  successful  in  the  acoustic  scenario.  For  example,  using  a 
linear  array  of  acoustic  sources,  and  moving  the  array  across  the  soil  surface,  it  should  be  possible 
to  obtain  high  lateral  resolution  in  the  direction  of  travel  of  the  array.  The  fractional  SAR  imaging 
algorithm  must  be  modified,  however,  to  account  for  signal  attenuation  in  the  soil,  and  to  counter 
the  effects  of  noise.  For  near-field  imaging,  we  plan  to  investigate  a  modified  SAR  imaging 
algorithm. 
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In  the  acoustic  subsurface  imaging  scenario,  the  “object”  to  be  imaged  (i.e.,  cultural 
artifacts  in  soil)  will  literally  abut  the  transducer,  probably  through  a  water  interface.  Thus, 
imaging  at  small  depths  will  take  place  in  the  near  field  of  the  transducer.  In  the  case  of  near-field 
imaging,  the  Fourier  approximation  breaks  down,  so  that  the  typical  correlation-based  imaging 
algorithm  for  SAR  will  fail.  For  near-field  SAR-based  imaging  we  suggest  a  modified  version  of 
the  so-called  wavenumber  approach  (Cafforio  et  al.,  1991;  Choi  and  Munson,  submitted).  The 
wavenumber  algorithm  is  equivalent  to  the  wave  migration  method  in  seismic  data  processing. 

This  technique  has  the  advantage  that  wavefront  curvature  is  incorporated  into  the  imaging  model 
so  that  superb  image  reconstructions  can  be  produced  at  close  ranges.  The  computational  demands 
of  this  algorithm  are  comparable  to  those  for  the  correlation-based  approach. 

The  SAR  technique  may  be  confounded  by  noise  introduced  in  the  subsurface  imaging 
scenario.  To  combat  this  potential  difficulty  we  can  utilize  a  new  SAR  imaging  algorithm  (Lee  et 
al.,  1996)  derived  from  a  geophysics  approach  used  to  estimate  the  conductivity  distribution  from 
wellbore  induction  measurements.  This  approach  to  SAR  imaging  implements  a  full  nonseparable 
two-dimensional  inversion  of  the  SAR  measurement  equation.  Using  the  spectral  representation 
of  the  SAR  measurement  kernel,  this  approach  decomposes  the  SAR  imaging  problem  into  a  set  of 
integral  equations,  described  by  a  nonseparable  convolution  in  one  dimension  and  a  projection  in 
the  other  dimension.  The  inversion  is  performed  in  two  stages:  multichannel  deconvolution, 
followed  by  back-projection.  By  dropping  small  singular  values  in  the  spectral  representation  and 
also  using  a  regularized  deconvolution  filter,  this  approach  to  SAR  imaging  is  inherently  robust  to 
the  presence  of  noise. 

A  third  image  processing  strategy  is  the  same  conventional  B-scan  imaging.  The  synthetic 
aperture  processing  offers  the  advantage  that  the  summed  beams  can  be  focused  at  any  depth; 
therefore,  the  entire  image  can  be  in  focus.  However,  the  disadvantage  of  synthetic  aperture 
imaging  is  that  the  signal  to  noise  ratio  is  low.  These  two  techniques  will  be  combined  and 
evaluated.  At  lower  acoustic  frequencies,  with  a  wider  beam  pattern  and  a  correspondingly  smaller 
aperture,  a  synthetic  approach  is  more  likely  whereas  at  higher  acoustic  frequencies,  with  a  more 
collimated  beam  pattern,  a  B-Scan  approach  is  more  logicd. 

2.  Experimental  component 

Professor  O’Brien  attended  the  Ultrasound  Transducer  Workshop  at  Penn  State,  August 
16-18, 1995,  as  an  invited  speaker.  There  he  made  valuable  contacts  with  engineers  and  scientists 
at  Penn  State’s  Applied  Research  Laboratory.  This  laboratory  has  a  major  activity  in  the  design 
and  development  of  high-power  sonar  transducers  in  our  acoustic  frequency  range,  mostly  arrays. 
The  cooperation  with  Penn  State’s  Applied  Research  Laboratory  provides  valuable  expertise  with 
the  two-dimensional  acoustic  transduction  device  and  some  of  the  support  electronics  including  RF 
signal  access  to  the  electrical  signals  received  at  each  of  the  52  elements.  In  addition,  significant 
expertise  exists  in  the  quantitative  and  qualitative  evaluation  of  the  imaging  system’s  overall 
development  and  performance  (see  Figure  18). 

A  52-element  sonar  array  (8x8  3.56  cm^  close-packed  elements  with  3  elements  in  each 
comer  missing)  can  be  driven  in  phase  to  produce  the  transmit  acoustic  source.  This  procedure 
will  produce  essentially  a  transmit  plane  wave.  This  is  a  cost-effective  means  to  evaluate  the 
feasibility  of  subsurface  imaging.  A  more  complete  (and  costly)  system  could  include  individual 
control  of  each  transmit  element,  thus  allowing  for  transmit  focusing  and  beam  steering. 
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On  receive,  the  backscattered  echo  signals  from  each  of  the  52  elements  will  be  amplified 
and  digitized.  Transmit  electronics  are  available  from  a  previously  funded  CERL  Project 
(DACA88-94-D-0008).  The  sonar  array  with  transmit^eceive  switches  and  receive  preamps  on 
each  element  can  be  modified  by  Penn  State’s  Applied  Research  Laboratoiy.  The  r^eive  digitizing 
capability  will  be  purchased  as  part  of  this  project  and  consists  of  one  Data  Translation  DT2839 
A/D  boards  capable  of  100  MHz  throughput  with  the  existing  133  MHz  Pentium  PC  in 
combination  with  a  DT  2896  Channel  Expander  capable  of  switching  betw^n  and  digitizing  the 
received  echoes  from  the  52  elements.  Tttis  system  operates  at  a  12-bit  digitizing  rate  of  1  MHz  to 
provide  the  operational  A/D  rate  per  channel  up  to  1  MHz.  Therefore,  this  cost-effective  A/D 
arrangement  will  collect  backscattered  echo  signals  from  all  52  elements.  This  proposed  approach 
will  require  multiple  (52)  transmit-receive  pulses  since  the  digitizing  capability  wiU  acquire  one 
echo  for  each  transmit  pulse.  Again,  this  is  a  cost-effective  means  to  evaluate  the  imaging 
feasibility.  A  more  complete  system  would  include  the  capability  to  a^uirc  the  electrical  signals 
from  each  transducer  element  for  a  single  transmit  pulse,  tiius  decreasing  the  imaging  data 
acquisition  time  and  providing  the  capability  to  average  received  signals  to  improve  the  signal-to- 
noise  ratio. 

Digitizing  the  receive  echo  signals  permits  off-line  processing  for  dynamic  focusing  on 
receive.  Table  5  provides  one-way  receive  focusing  properties  of  the  52-element  sonar  array 
assuming  an  active  aperture  diameter  of  28  cm  in  a  medium  with  a  propagation  speed  of  250  m/s. 
The  regions  denoted  by  dashes  in  Table  5  do  not  permit  focusing  under  the  indicated  conditions 
(focusing  must  occur  in  the  Fresnel  region  of  the  field)  and  thus,  in  these  regions,  synthetic 
aperture  processing  is  more  logical  to  optimize  lateral  resolution. 


Table  5.  Beam  properties  for  a  28  cm  active  aperture  diameter  in  a  medium  with  a  propagation 
speed  of  250  m/s.  A  20%  bandwidth  is  assumed  for  axial  resolution. 


Aunroximate  Lateral  Resolution  at  Focus  (cm) 
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5.0 

— 

— 

— 

0.25  1 
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55 
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— 

— 

0.18 

8 

28 
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— 

0.16 
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32 
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4.2 
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10 

35 
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78 

2.5 

3.8 

5.0 
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11 

39 

2.3 

86 

2.3 

3.4 

4.5 

5.7 

0.11 

12 

42 

2.1 

94 

2.1 

3.1 

4.2 

5.2 

0.10 

13 

46 

1.9 
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1.9 

2.9 

3.8 

4.8 

0.096 

14 

49 

1.8 

110 

1.8 

2.7 

3.6 

0.089 

15 

53 

1.7 

118 

1.7 

2.5 

3.3 

4.2 

0.083 

At  a  digitizing  rate  of  200  kHz  (at  least  10  times  greater  that  the  highest  imaging 
frequency),  and  assuming  a  soil  propagation  speed  of  250  m/s,  400  bytes  of  data  are  collected  per 
meter  depth  (denoted  a  single  A-line).  Assuming  1  kB/A-line  (2.5  m  depth;  approximately  8.2 
feet),  data  acquisition  from  one  array  position  yields  52  digitiz^  A-lines  and  52  kB/position. 
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The  airay  can  be  moved  with  two  geometries,  one  in  a  straight  line  to  yield  a  2-D  image 
(denoted  B-mode  image)  and  the  other  in  a  2-D  pattern  to  yield  a  3-D  image.  The  spacing  between 
array  positions  should  be  about  2  cm  to  take  full  advantage  of  independent  A-lines  (note:  beam 
diameter  at  focus  is  about  2  cm  at  15  kHz  for  f/1  receive  focusing,  see  Table  5)  for  B-mode 
image  processing  and  possibly  a  higher  density  (<2mm  spacing)  for  synthetic  array  image 
processing.  About  10.4  MB  of  data  (200x52  kB)  are  collected  for  2-D,  B-mode  imaging  with  200 
array  positions  (potential  2-D  field  of  view:  4  m  horizontally  x  2.5  m  in  depth),  whereas  about  520 
MB  of  data  (100x100x52  kB)  are  collected  for  3-D  imaging  with  a  100  x  100  square  2-D  pattern 
(potential  3-D  field  of  view:  2  x  2  m  horizontally  x  2.5  m  in  depth).  The  advantage  of  the  3-D 
imaging  data  acquisition  procedure  is  that  not  only  multiple  2-D  images  perpendicular  to  the  soil 
surface  can  be  generated  but  also  multiple  2-D  horizontal  images  parcel  to  the  soil  surface 
(denoted  C-mcde  images)  at  any  soil  depth  can  be  generated. 

Test  bed  recommendations  It  is  recommended  that  CERL’s  test  imaging  field  include  a 
series  of  2  cm  and  5  cm  diameter  metal  pipes  separated  at  various  distances  horizontally  and 
vertically  and  at  various  depths  to  provide  high  contrast  targets  so  tiiat  axial  and  lateral  resolutions 
can  be  estimated  experimentally.  Further,  it  is  recommended  that  CERL’s  test  imaging  field 
include  a  similar  series  of  thin-walled  plastic  containers  filled  with  materials  of  v^ous  acoustic 
scattering  properties  to  estimate  experimentally  contrast  resolution.  A  3-D  positioning  system  will 
be  required  to  support  the  300-pound  array  and  be  capable  of  moving  the  array  in  about  2  cm 
increments  with  a  precision  of  about  2  mm  (about  a  factor  of  10  better  than  the  2  cm  incremental 
movement).  The  top  surface  of  the  test  bed  should  be  water  (depth  about  6-8  inches)  to  provide 
acoustic  coupling  between  the  array  source  and  soil.  Finally,  the  array  positioning  system  should 
be  interfaced  with  a  computer  to  precisely  move  and  record  the  position.  If  possible,  it  would  be 
desirable  to  have  the  positioning  system  interfaced  to  the  data  acquisition  system.  Also,  the  drive 
electronics  for  the  52-element  array  require  a  208  volt,  three-phase  line. 

E.  FUTURE  PLANS 

The  long-term  objective  of  the  proposed  research  program  is  to  classify  precisely 
subsurface  cultural  artifacts  through  acoustic  imaging  means.  Six  generalized  research  phases  are 
envisioned  to  accomplish  the  long-term  objective.  Phase  1  is  the  currently  funded  activity 
pro^am  (CERL  Contract  Number  DACA88-94-D-0008,  “Acoustical  Characterization  of  Soil”) 
and  is  providing  base-line  evaluations  of  acoustic  properties  of  a  limited  set  of  soils.  Phase  2 
evaluates  the  feasibility  for  subsurface  detection  of  cultural  artifacts  using  various  acoustic  imaging 
approaches  under  laboratory,  test-bed  conditions  and  to  initiate  studies  of  in  situ  soil  acoustic 
properties.  Phase  3  optimizes  both  the  imaging  system  and  the  image  formation  algorithms  for 
detection  of  cultural  artifacts  under  laboratory,  test-bed  conditions.  Phase  4  develops 
classification  strategies  to  classify  precisely  subsurface  cultural  artifacts  under  laboratory,  test-bed 
conditions,  and  to  test  detection  strategies  under  in  situ,  field  conditions.  Phase  5  evaluates  and 
optimizes  detection  and  classification  strategies  under  selected  in  situ,  field  conditions  This  phase 
consists  of  developing  the  acoustic  imaging  system  on  an  appropriate  mobile  platform.  Phase  6 
evaluates  detection  and  classification  capabilities  under  a  wide  variety  of  in  situ,  field  conditions, 
and  modifies  these  capabilities  as  appropriate. 

All  of  these  phases  require  continued  evaluation  of  the  acoustic  propagation  and  scattering 
properties  of  soils.  Thus,  it  is  proposed  to  continue  these  fundamental  acoustic  characterizations  of 
well-characterized  soil  samples  to  develop  a  better  understanding  of  the  acoustical  properties  under 
more  typical  soil  conditions.  The  acoustic  imaging  feasibility  evaluation  in  CERL’s  test  bed  will  be 
conducted  in  one  type  of  soil  whereas  there  are  numerous  and  widely  varying  soil  conditions 
around  the  world.  The  increased  understanding  of  soil’s  acoustic  properties  will  permit  a 
fundamental  basis  for  extrapolating  imaging  capabilities  in  soil  types  that  have  not  been  evaluated 
in  a  soil  test  bed. 
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It  is  proposed  to  add  carbonate  content,  plastic  limit,  liquid  limit,  maximum  density  and 
coefficient  of  linear  extensibility  (COLE)  testing  to  the  analyses  currently  being  done  on  the  soil 
samples  currently  under  evaluation.  Further,  it  is  proposed  to  drop  the  soil  characteristics  that  do 
not  correlate  with  acoustic  properties  and  re-evaluate  using  those  that  do  correlate  along  with  the 
additional  soil  characteristics.  These  analyses  are  critical  for  assessing  whether  the  soil  properties 
can  provide  a  fundamental  basis  for  estimating  the  soil’s  acoustic  properties  and  hence  the  acoustic 
imaging  capabilities.  The  additional  analyses  will  be  conducted  on  the  current  6  soil  types. 

A  second  component  to  be  investigated  is  the  effect  of  soil  layering  which  can  be  performed 
with  the  soil  samples  that  have  already  been  collected  and  evaluated,  and  with  the  existing 
measurement  system.  Here,  the  layering  would  be  accomplished  with  a  thin  plastic  film  separating 
the  layers.  The  soil  layers  in  a  single  study  would  include  different  soils  (or  the  same  soil  with 
different  water  contents  or  compactions)  where  the  selection  of  the  specific  soils  would  be  based 
on  their  individual  acoustic  properties.  Mathematical  modeling  would  be  employed  to  determine 
whether  the  acoustic  properties  of  individual  soil  samples  can  estimate  the  acoustic  properties  of  the 
layered  media. 

A  third  component  to  be  investigated  is  undisturbed  core  samples.  A  coring  truck  exists  to 
collect  about  6”  diameter  core  samples  about  24”  deep.  This  is  an  ideal  size  for  the  existing 
measurement  system.  The  intent  would  be  to  evaluate  undisturbed  core  samples  of  soil  types  that 
have  already  been  evaluated  and  thus  to  compare  the  differences  between  the  same  soil  typ;s  that 
have  been  prepared  under  laboratory  conditions  and  those  that  are  more  typical  of  undisturbed  in 
situ  conditions.  Ten  paired  core  samples  will  be  obtained,  one  for  acoustic  characterization  and  the 
other  for  soil  characterization. 
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Appendix  Al.  Soil  conditions  at  time  of  acoustic  evaluation. 
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2.0 

16.0 

17.10 

0.80 

0.92 

2.49 

32.2 

67.8 

17.0 

11.5 

14.3 

12.5 

adal9211-D2 

20.2 

17.2 

3.0 

19.0 

20.43 

0.84 

0.99 

2.49 

33.8 

66.2 

22.1 

14.7 

17.4 

14.8 

ada22211-Dl 

23.6 

20.1 

3.5 

22.0 

23.76 

0.85 

0.99 

2.49 

33.9 

66.1 

22.3 

14.7 

17.4 

14.8 

adal0241-El 

11.7 

9.9 

1.7 

10.0 

10.44 

0.95 

1.12 

2.49 

38.1 

61.9 

26.8 

16.5 

17.4 

14.8 

adal3241-E2 

16.1 

14.1 

2.0 

13.0 

13.77 

1.02 

1.17 

2.49 

41.0 

59.0 

24.8 

14.6 

14.3 

12.5 

adal6241-E3 

19.5 

17.1 

warn 

16.0 

17.10 

1.00 

1.14 

2.49 

40.0 

60.0 

23.8 

14.3 

14.3 

12.5 

adal9241-E4 

23.7 

20.7 

3.0 

19.0 

20.43 

1.01 

1.16 

2.49 

40.6 

59.4 

24.4 

14.5 

14.3 

12.5 

adal6311-F4 

17.2 

13.2 

4.0 

16.0 

17.10 

0.77 

1.01 

2.49 

30.9 

69.1 

34.2 

23.7 

30.8 

23.5 

adal9311-F3 

21.1 

16.1 

5.0 

19.0 

20.43, 

0.79 

1.03 

2.49 

31.7 

68.3 

35.6 

24.3 

30.8 

23.5 

ada22311-F2 

24.8 

19.0 

5.8  . 

22.0 

23.76 

0.80 

1.04 

2.49 

32.0 

68.0 

36.1 

24.6 

30.8 

23.5 

ada25311-Fl 

28.0 

21.6 

25.0 

27.09 

0.80 

1.03 

2.49 

32.0 

68.0 

34.9 

23.7 

29.8 

23.0 

adal0341-G4 

12.3 

EEI 

2.9 

10.0 

10.48 

0.90 

1.17 

2.49 

36.0 

64.0 

43.2 

27.7 

30.8 

EBB 

adal3341-G3 

16.3 

12.5 

3.8 

13.0 

13.81 

0.90 

1.18 

2.49 

36.2 

63.8 

43.5 

27.8 

30.8 

23.5 

adal6341-G2 

21.1 

16.1 

5.0 

16.0 

17.14 

0.94 

1.23 

2.49 

37.8 

62.2 

46.6 

29.0 

30.8 

23.5 

adal9341-Gl 

24.6 

18.8 

5.8 

19.0 

20.47 

0.92 

1.20 

2.49 

36.9 

63.1 

44.8 

28.3 

30.8 

23.5 

adal6511-H3 

24.8 

14.7 

10.1 

15.0 

15.99 

0.92 

1.55 

2.49 

36.8 

63.2 

100.1 

63.2 

68.9 

40.8 

adal9511-H4 

26.5 

15.7 

10.8 

16.0 

17.10 

0.92 

1.55 

2.49 

36.8 

63.2 

100.1 

63.2 

68.9 

40.8 

ada20511-H6 

33.4 

19.8 

13.6 

20.0 

21.54 

0.92 

1.55 

2.49 

36.8 

63.2 

ngiin 

63.2 

68.9 

40.8 

ada22511-H2 

36.8 

21.8 

15.0 

22.0 

23.76 

0.92 

1.55 

2.49 

36.8 

63.2 

100.1 

63.2 

68.9 

40.8 

ada23511-H5 

38.5 

22.8 

15.7 

23.0 

24.87 

0.92 

1.55 

2.49 

36.8 

63.2 

100.1 

63.2 

68.9 

40.8 

ada25511-Hl 

42.0 

24.9 

17.1 

25.0 

27.09 

0.92 

1.55 

2.49 

36.8 

63.2 

100.1 

63.2 

68.9 

40.8 
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Appendix  Al.  Soil  conditions  at  time  of  acoustic  evaluation  (continued). 


Sample  File 
Number 

Net 

Total 

Dry 

Soil 

kp 

Net 

HP 

5^8 

Tub 

Depth 

cm 

Total 

Vol 

L 

Dry 

Soil 

Den 

gm/cc 

Total 

Wet 

Den 

gm/cc 

Soil 

Part 

Den 

Rm/cc 

Total 

Solid 

% 

Total 

Pores 

% 

HP 

Filled 

Pores 

% 

% 

HP 

(vol) 

% 

HP 

(dry) 

% 

Hp 

(wet) 

cabl6111-A4 

23.6 

22.8 

0.8 

16.0 

16.97 

1.34 

1.39 

2.64 

50.9 

49.1 

mm 

3.3 

cabl9111-A3 

28.7 

27.7 

1.0 

19.0 

20.30 

1.37 

1.41 

2.64 

51.7 

48.3 

msM 

m 

wm 

3.3 

cab22111-A2 

32.7 

31.6 

1.1 

22.0 

23.63 

1.34 

1.38 

2.64 

50.5 

49.5 

9.3 

E&l 

warn 

3.3 

cab25111-Al 

38.6 

37.3 

1.3 

25.0 

26.96 

1.38 

1.43 

2.64 

52.3 

47.7 

10.0 

EO 

mm 

3.3 

cabl6112-Bl 

19.4 

18.7 

0.6 

16.0 

17.17 

1.09 

1.13 

2.64 

41.2 

58.8 

mm 

Bl 

3.3 

cabl9112-B2 

24.9 

24.1 

0.8 

19.0 

20.50 

1.17 

1.21 

2.64 

44.4 

55.6 

4.0 

Rl 

3.3 

cab22112-B3 

29.7 

28.7 

1.0 

22.0 

23.83 

1.20 

1.25 

2.64 

45.6 

54.4 

mm 

4.1 

El 

3.3 

cab25112-B4 

34.7 

33.5 

1.2 

25.0 

27.16 

1.24 

1.28 

2.64 

46.7 

53.3 

8.0 

Ea 

El 

3.3 

cabl0141-Cl 

IS.O 

14.5 

0.5 

10.0 

10.31 

1.41 

1.45 

2.64 

53.2 

46.8 

10.3 

El 

3.3 

cabl3141-C2 

20.0 

19.3 

0.7 

13.0 

13.64 

1.41 

1.46 

2.64 

53.5 

46.5 

10.4 

EO 

El 

3.3 

cabl6141-C3 

2S.7 

24.8 

0.9 

16.0 

16.97 

1.46 

1.51 

2.64 

55.4 

44.6 

11.3 

5.0 

El 

3.3 

cabl9141-C4 

30.6 

29.6 

1.0 

19.0 

20.30 

1.46 

1.51 

2.64 

55.1 

44.9 

11.2 

5.0 

El 

3.3 

cabl6211-D4 

21.6 

19.8 

1.8 

16.0 

16.97 

1.17 

1.27 

2.64 

44.2 

55.8 

18.5 

10.3 

8.8 

8.1 

cabl9211-D3 

25.3 

23.2 

2.1 

19.0 

20.30 

1.15 

1.25 

2.64 

43.3 

56.7 

17.9 

10.1 

8.8 

8.1 

cab22211-D2 

30.0 

27.7 

2.3 

22.0 

23.63 

1.17 

1.27 

2.64 

44.3 

55.7 

17.5 

9.8 

8.3 

WXM 

cab25211-Dl 

33.9 

31.3 

2.6 

25.0 

26.96 

1.16 

1.26 

2.64 

43.9 

56.1 

17.2 

wm 

8.3 

El 

cabl0241-El 

14.4 

13.2 

1.2 

10.0 

10.27 

1.29 

1.40 

2.64 

48.7 

51.3 

22.2 

11.4 

8.8 

8.1 

cabl3241-E2 

19.2 

17.6 

1.6 

13.0 

13.57 

1.30 

1.41 

2.64 

49.2 

50.8 

22.6 

11.5 

8.8 

8.1 

cabl6241-E3 

25.7 

23.6 

2.1 

16.0 

16.87 

1.40 

1.52 

2.64 

53.0 

47.0 

26.3 

12.4 

8.8 

8.1 

cabl9241-E4 

31.9 

29.3 

2.6 

19.0 

20.17 

1.45 

1.58 

2.64 

55.0 

45.0 

28.5 

12.8 

8.8 

8.1 

cabl6311-F4 

17.2 

14.5 

WSM 

16.0 

17.01 

0.85 

1.01 

2.64 

32.2 

67.8 

23.4 

15.9 

18.6 

15.7 

cabl9311-F3 

21.2 

17.9 

3.3  ■ 

19.0 

20.34 

0.88 

1.04 

2.64 

33.4 

66.6 

24.0 

16.0 

18.1 

15.4 

cab22311-F2 

24.3 

20.6 

EEi 

22.0 

23.67 

0.87 

1.03 

2.64 

32.9 

67.1 

23.5 

15.8 

18.1 

15.4 

cab25311-Fl 

29.1 

24.6 

Ea 

25.0 

27.00 

0.91 

1.08 

2.64 

65.5 

25.3 

16.5 

18.1 

15.4 

cabl0341-G4 

12.7 

10.7 

2.0 

10.0 

10.35 

1.03 

1.23 

2.64 

39.1 

60.9 

31.7 

19.3 

18.6 

15.7 

cabl3341-G3 

19.0 

16.0 

3.0 

13.0 

13.68 

1.17 

1.39 

2.64 

44.3 

55.7 

39.2 

21.8 

18.6 

15.7 

cabl6341-G2 

22.9 

19.3 

3.6 

16.0 

17.01 

1.13 

1.35 

2.64 

42.9 

57.1 

37.1 

21.2 

18.6 

15.7 

cabl9341-Gl 

29.6 

24.9 

in 

19.0 

20.34 

1.23- 

1.46 

2.64 

46.4 

53.6 

42.7 

22.9 

18.6 

15.7 

cables  11-H6 

17.7 

11.6 

6.1 

10.0 

10.46 

1.11 

1.69 

2.64 

42.0 

58.0 

100.0 

58.0 

52.2 

34.3 

cabl3511-H5 

23.3 

15.3 

8.0 

13.0 

13.79 

1.11 

1.69 

2.64 

42.0 

58.0 

100.0 

58.0 

52.2 

cables  11-H4 

28.9 

19.0 

9.9 

16.0 

17.12 

1.11 

1.69 

2.64 

42.0 

58.0 

100.0 

58.0 

52.2 

34.3 

cabl7Sll-H8 

30.8 

20.2 

10.6 

17.0 

18.23 

1.11 

1.69 

2.64 

42.0 

58.0 

100.0 

58.0 

52.2 

34.3 

cabl9Sll-H3 

34.6 

22.7 

11.9 

19.0 

20.45 

1.11 

1.69 

2.64 

42.0 

58.0 

100.0 

58.0 

52.2 

34.3 

cab21Sll-H7 

38.3 

25.2 

13.1 

21.0 

22.67 

1.11 

1.69 

2.64 

42.0 

58.0 

100.0 

58.0 

52.2 

34.3 

cab22Sll-H2 

40.2 

26.4 

22.0 

23.78 

1.11 

1.69 

2.64 

42.0 

58.0 

100.0 

58.0 

52.2 

34.3 

cab2SSll-Hl 

45.8 

15.7 

25.0 

27.11 

1.11 

1.69 

2.64 

42.0 

58.0 

100.0 

58.0 

52.2 

34.3 
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Appendix  Al.  Soil  conditions  at  time  of  acoustic  evaluation  (continued). 


Sample  File 
Number 


Total  S( 
Wet  P 
Den 
em/cc 


Total 

Total 

HP 

Filled 

Solid 

Pores 

Pores 

% 

% 

% 

dral2111-A4 

17.0 

16.2 

0.7 

12.0 

12.46 

1.30 

1.36 

2.52 

51.8 

48.2 

11.9 

EEIKEIESI 

dral7112-A3 

23.9 

22.8 

1.0 

17.0 

17.91 

1.28 

1.33 

2.52 

50.7 

49.3 

11.4 

5.6 

mxM 

mm 

dra22113-A2 

31.6 

30.3 

1.3 

22.0 

23.36 

1.30 

1.35 

2.52 

51.5 

48.5 

11.8 

En 

mm 

dra27112-Al 

40.2 

38.5 

1.7 

27.0 

28.81 

1.33 

1.39 

2.52 

53.0 

47.0 

12.5 

E9 

EEI 

mm 

dral2114-B4 

18.8 

18.0 

0.8 

12.0 

12.75 

1.41 

1.47 

2.52 

56.1 

43.9 

14.2 

6.2 

warn 

Ea 

dral7114-B3 

25.2 

24.1 

1.1 

17.0 

lEBH 

1.32 

1.37 

2.52 

52.3 

47.7 

12.2 

5.8 

wm 

mm 

dra22114-B2 

31.7 

30.3 

1.3 

22.0 

23.85 

1.27 

1.33 

2.52 

50.5 

49.5 

11.3 

5.6 

mm 

mm 

dra27114-Bl 

40.6 

38.8 

1.7 

27.0 

29.40 

1.32 

1.38 

2.52 

52.5 

47.5 

12.3 

5.8 

mm 

Ka 

dral2141-Cl 

17.0 

16.2 

0.7 

12.0 

12.46 

1.30 

1.36 

2.52 

51.8 

48.2 

11.9 

mm 

mm 

dral7141-C2a 

27.8 

26.7 

1.2 

17.0 

17.91 

1.49 

1.55 

2.52 

59.1 

40.9 

16.1 

6.6 

KEI 

Ka 

dral7142-ab 

27.8 

26.7 

1.2 

17.0 

17.91 

1.49 

1.55 

2.52 

59.1 

40.9 

16.1 

6.6 

mm 

mm 

dra23141-C3 

37.0 

35.4 

1.6 

23.0 

24.45 

1.45 

1.51 

2.52 

57.5 

42.5 

15.0 

mm 

mm 

dra05143-Dl 

6.3 

6.0 

0.3 

5.0 

4.83 

1.24 

1.29 

2.52 

49.2 

50.8 

10.8 

5.5 

mm 

mm 

dra08143-D2 

12.0 

11.4 

0.5 

8.0 

8.10 

1.41 

1.48 

2.52 

56.1 

43.9 

14.2 

6.2 

KEI 

mm 

dral0143-D3 

15.8 

15.1 

0.7 

10.0 

10.28 

1.47 

1.53 

2.52 

58.3 

41.7 

15.5 

6.5 

KEIEU 

dral6211-E4 

20.3 

17.7 

2.6 

16.0 

16.88 

1.05 

1.20 

2.52 

41.6 

58.4 

26.1 

15.3 

14.6 

M.l 

dral9211-E3 

24.1 

21.0 

3.1 

19.0 

20.15 

1.04 

1.19 

2.52 

41.4 

58.6 

25.9 

15.2 

14.6 

12.7 

dra22211-E2 

28.4 

24.8 

E9 

22.0 

23.42 

1.06 

1.21 

2.52 

42.0 

58.0 

26.6 

15.4 

14.6 

12.7 

dra25211-El 

31.7 

27.7 

4.0 

25.0 

26.69 

1.04 

1.19 

2.52 

41.2 

58.8 

25.7 

15.1 

14.6 

12.7 

dral6212-Fl 

17.9 

15.6 

2.3 

16.0 

16.91 

0.92 

1.06 

2.52 

36.7 

63.3 

21.3 

13.5 

14.6 

12.7 

dral9212-F2 

21.3 

18.6 

WOk 

19.0 

20.21 

0.92 

1.05 

2.52 

36.5 

63.5 

21.1 

13.4 

14.6 

12.7 

dra22212-F3 

25.7 

22.4 

3.3 

22.0 

23.51 

0.95 

1.09 

2.52 

37.9 

62.1 

22.4 

13.9 

14.6 

12.7 

dra25212-F4 

30.2 

26.4 

3.8 

25.0 

26.81 

0.98 

1.13 

2.52 

39.1 

60.9 

23.5 

14.3 

14.6 

12.7 

dral0241-Gl 

13.8 

12.1 

1.7 

10.0 

10.35 

1.17 

1.33 

2.52 

46.6 

53.4 

30.1 

16.1 

13.7 

12.1 

dral3241-G2 

18.5 

16.3 

2.2 

13.0 

13.68 

1.19 

1.35 

2.52 

47.2 

52.8 

30.9 

16.3 

13.7 

12.1 

dral6241-G3 

22.8 

20.1 

16.0 

17.01 

1.18 

1.34 

2.52 

46.8 

53.2 

30.4 

16.2 

13.7 

12.1 

dral9241-G4 

28.5 

25.1 

Bl 

19.0 

20.34 

1.23 

1.40 

2.52 

49.0 

51.0 

33.1 

16.9 

13.7 

12.1 

dral6311-H4 

18.2 

14.4 

3.8 

16.0 

16.91 

0.85 

1.07 

2.52 

33.8 

66.2 

33.8 

22.4 

20.8 

dral9311-H3 

23.3 

18.0 

5.3 

19.0 

20.21 

0.89 

1.15 

2.52 

35.4 

64.6 

40.4 

26.1 

29.3 

22.6 

dra22311-H2 

27.0 

20.9 

6.1 

22.0 

23.51 

0.89 

1.15 

2.52 

35.2 

64.8 

40.1 

25.9 

29.3 

22.6 

dra25311-Hl 

31.1 

24.1 

7.0 

25.0 

26.81 

0.90 

1.16 

2.52 

35.7 

64.3 

40.8 

26.3 

29.3 

22.6 

dral0341-Il 

15.1 

12.0 

3.1 

10.0 

10.31 

1.16 

1.46 

2.52 

46.1 

53.9 

56.6 

30.5 

26.3 

20.8 

dral3341-I2 

22.5 

17.8 

Kn 

13.0 

13.61 

1.31 

1.65 

2.52 

51.9 

48.1 

71.4 

34.4 

26.3 

20.8 

dral6341-B 

26.9 

21.2 

msM 

16.0 

16.91 

1.25 

1.59 

2.52 

49.7 

50.3 

67.5 

33.9 

27.1 

21.3 

dral8341-I4 

30.7 

24.2 

6.5 

17.5 

18.56 

1.30 

1.65 

2.52 

51.7 

48.3 

73.0 

35.3 

27.1 

21.3 

dral5511-J8 

26.9 

18.0 

8.8 

15.0 

15.99 

1.13 

1.68 

2.52 

44.8 

55.2 

100.1 

55.3 

49.0 

32.9 

dral6511-J4 

28.7 

19.3 

9.5 

16.0 

17.10 

1.13 

1.68 

2.52 

44.8 

55.2 

100.1 

55.3 

49.0 

32.9 

dral8511-J7 

32.5 

21.8 

10.7 

18.0 

19.32 

1.13 

1.68 

2.52 

44.8 

55.2 

100.1 

55.3 

49.0 

32.9 

dral9511-J3 

34.3 

23.0 

11.3 

19.0 

20.43 

1.13 

1.68 

2.52 

44.8 

55.2 

100.1 

55.3 

49.0 

32.9 

dra21511-J6 

38.1 

25.5 

12.5 

21.0 

22.65 

1.13 

1.68 

2.52 

44.8 

55.2 

100.1 

55.3 

49.0 

32.9 

dia22511-J2 

39.9 

26.8 

13.1 

22.0 

23.76 

1.13 

1.68 

2.52 

44.8 

55.2 

100.1 

55.3 

49.0 

32.9 

dra24511-J5 

43.6 

29.3 

14.4 

24.0 

25.98 

1.13 

1.68 

2.52 

44.8 

55.2 

100.1 

55.3 

49.0 

32.9 

dra25511-Jl 

45.5 

30.5 

15.0 

25.0 

27.09 

1.13 

1.68 

2.52 

44.8 

55.2 

100.1 

55.3 

49.0 

32.9 
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Appendix  Al.  Soil  conditions  at  time  of  acoustic  evaluation  (continued). 


Sample  File 
Number 


meal6111-A4 


meal9111-A3 


mea22111-A2 


mea25111-Al 


meal6112-Bl 


meal9112-B2 


mea22112-B3 


mea25112-B4 


meal0141-Cl 


meal3141-C2 


meal6141-C3 


meal9141-C4 


meal6211-D4 


meal9211-D3 


mea22211-D2 


mea25211-Dl 


meal0241-El 


meal  324 1-E2 


meal6241.E3 


meal9241-E4 


meal6311-F4 


meal9311-F3 


mea22311-F2 


mea25311-Fl 


meal0341-Gl 


meal3341-G2 


meal6341-G3 


meal9341-G4 


meal0511-H6 


meal2511-H5 


eal3511-H10 


meal5511-H4 


meal6511-H9 


meal8511-H3 


meal9511-H8 


mea21511-H2 


mea22511-H7 


mea24511-Hl 


Tub 

Depih 

cm 


16.0 


19.0 


22.0 


25.0 


16.0 


19.0 


22.0 


25.0 


10.0 


13.0 


16.0 


19.0 


16.0 


19.0 


22.0 


25.0 


10.0 


13.0 


16.0 


19.0 


16.0 


19.0 


22.0 


25.0 


10.0 


13.0 


16.0 


19.0 


10.0 


12.0 


13.0 


15.0 


16.0 


18.0 


19.0 


21.0 


22.0 


24.0 


Total 

Total 

HP 

Filled 

Solid 

Pores 

Pores 

% 

% 

% 

17.12 


20.45 


23.78 


27.11 


17.14 


20.47 


23.80 


27.13 


10.46 


13.79 


17.12 


20.45 


17.12 


20.45 


23.78 


27.11 


10.31 


13.61 


16.91 


20.21 


16.87 


20.17 


23.47 


26.77 


10.36 


13.69 


17.02 


20.35 


10.31 


12.51 


13.61 


15.81 


16.91 


19.11 


20.21 


22.41 


23.51 


25.71 
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Appendix  Al.  Soil  conditions  at  time  of  acoustic  evaluation  (continued). 


Sample  File 
Number 


J 

)lalllll-A6 

17.7 

17.6 

0.1 

11.0 

11.37 

1.55 

1.56 

2.65 

58.7 

41.3 

1.1 

0.5 

0.3 

0.3 

J 

»lal4111-A5 

21.8 

21.7 

0.1 

14.0 

14.67 

1.48 

1.49 

2.65 

56.0 

44.0 

1.0 

0.4 

0.3 

0.3 

J 

)lal7111-A4 

26.3 

26.2 

0.1 

17.0 

17.97 

1.46 

1.46 

2.65 

55.2 

44.8 

0.9 

0.4 

0.3 

0.3 

J 

)la20111-A3 

32.6 

32.5 

0.1 

20.0 

21.27 

1.53 

1.53 

2.65 

57.8 

42.2 

1.1 

0.4 

0.3 

0.3 

J 

>la23111-A2 

37.1 

37.0 

0.1 

23.0 

24.57 

1.51 

1.51 

2.65 

56.9 

43.1 

1.0 

0.4 

0.3 

0.3 

J 

)la26111-Al 

41.2 

41.1 

0.1 

25.5 

27.32 

1.50 

1.51 

2.65 

56.8 

43.2 

1.0 

0.4 

0.3 

0.3 

J 

»lal2112-Bl 

17.8 

17.7 

0.1 

12.0 

12.67 

1.40 

1.40 

2.65 

53.0 

47.0 

0.9 

0.4 

0.3 

0.3 

J 

>lal7113-B2 

29.2 

29.1 

0.1 

17.0 

18.17 

1.60 

1.61 

2.65 

60.6 

39.4 

1.2 

0.5 

0.3 

0.3 

J 

)la22112-B3 

37.1 

37.0 

0.1 

22.0 

23.67 

1.56 

1.57 

2.65 

59.1 

40.9 

1.1 

0.5 

0.3 

0.3 

J 

Ia27112-B4 

48.1 

47.9 

0.1 

27.0 

29.17 

1.64 

1.65 

2.65 

62.1 

37.9 

1.3 

0.5 

0.3 

0.3 

J 

Ial2114-C4 

19.1 

19.0 

0.1 

12.0 

12.47 

1.53 

1.53 

2.65 

57.7 

42.3 

1.1 

0.4 

0.3 

0.3 

J 

Ial7114-C3 

27.1 

27.0 

0.1 

17.0 

18.28 

1.48 

1.48 

2.65 

55.9 

44.1 

1.0 

0.4 

0.3 

0.3 

J 

Ia22114-C2 

36.3 

36.1 

0.1 

22.0 

23.85 

1.52 

1.52 

2.65 

57.3 

42.7 

1.0 

0.4 

0.3 

0.3 

_E 

Ia27114-Cl 

44.5 

44.3 

0.1 

27.0 

29.38 

1.51 

1.51 

2.65 

57.0 

43.0 

1.0 

0.4 

0.3 

0.3 

J 

Ial2115-Dl 

17.5 

17.4 

0.1 

12.0 

12.73 

1.37 

1.37 

2.65 

51.7 

48.3 

0.8 

0.4 

0.3 

0.3 

J 

Ial7115-D2 

28.0 

27.9 

0.1 

17.0 

17.97 

1.55 

1.56 

2.65 

58.7 

41.3 

1.1 

0.5 

0.3 

0.3 

_E 

Ia22115-D3 

35.3 

35.2 

0.1 

22.0 

23.47 

1.50 

1.50 

2.65 

56.7 

43.3 

1.0 

0.4 

0.3 

0.3 

J 

Ia26115-D4 

43.9 

43.7 

0.1 

26.0 

27.87 

1.57 

BE3 

2.65 

59.3 

40.7 

1.1 

0.5 

0.3 

0.3 

_E 

Ia27115-D5 

45.6 

45.4 

0.1 

27.0 

28.97 

1.57 

1.57 

2.65 

59.3 

40.7 

1.1 

0.5 

0.3 

0.3 

_E 

Ial8216-E4 

31.6 

29.2 

EEI 

18.0 

19.07 

1.53 

1.66 

2.65 

57.9 

42.1 

29.8 

12.5 

8.2 

Mrici 

J 

Ia21216-E3 

37.9 

35.0 

2.9 

21.0 

22.37- 

1.57 

1.69 

2.65 

59.2 

40.8 

31.4 

12.8 

8.2 

BQ 

_E 

Ia24216-E2 

43.8 

40.5 

3.3  . 

24.0 

25.67 

1.58 

1.71 

2.65 

59.6 

40.4 

32.0 

12.9 

8.2 

_E 

la272 16-El 

50.1 

46.3 

3.8 

27.0 

28.97 

1.60 

1.73 

2.65 

60.4 

39.6 

33.0 

13.1 

8.2 

■Bi 

J 

Ial9217-Fl 

28.0 

25.4 

2.5 

19.0 

20.09 

1.26 

1.39 

2.65 

47.8 

52.2 

24.3 

12.7 

10.0 

9.1 

_E 

Ia22217-F2 

33.3 

30.2 

3.0 

22.0 

23.36 

1.29 

1.42 

2.65 

48.9 

51.1 

25.4 

13.0 

10.0 

9.1 

_E 

Ia25217-F3 

41.3 

37.5 

3.8 

25.0 

26.63 

1.41 

1.55 

2.65 

53.2 

46.8 

30.1 

14.1 

10.0 

9.1 

J 

Ial0241-Gl 

17.8 

16.4 

1.4 

10.0 

10.27 

1.60 

1.73 

2.65 

60.4 

39.6 

34.4 

13.6 

8.5 

_E 

Ial3241-G2 

23.2 

21.4 

1.8 

13.0 

13.57 

1.58 

1.71 

2.65 

59.6 

40.4 

33.2 

13.4 

8.5 

KEI 

_E 

Ial6241-G3 

29.9 

27.5 

16.0 

16.87 

1.63 

1.77 

2.65 

61.7 

38.3 

36.4 

13.9 

8.5 

_E 

Ial9241-G4 

33.6 

2.9 

19.0 

20.17 

1.67 

1.81 

2.65 

63.0 

37.0 

38.5 

14.2 

8.5 

KEI 

_E 

Ial6311-H4 

26.7 

23.9 

IBM 

16.0 

16.82 

1.42 

1.58 

2.65 

53.8 

46.2 

34.9 

16.2 

11.4 

10.2 

_E 

Ial9311-H3 

32.1 

29.6 

2.5 

19.0 

20.09 

1.47 

1.60 

2.65 

55.8 

44.2 

27.8 

12.3 

8.4 

_E 

Ia22311-H2 

37.3 

34.4 

2.9 

22.0 

23.36 

1.47 

1.59 

2.65 

55.6 

44.4 

27.7 

12.3 

8.4 

''WB 

_E 

Ia25311-Hl 

41.9 

38.6 

3.2 

25.0 

26.63 

1.45 

1.57 

2.65 

54.8 

45.2 

26.8 

12.1 

8.4 

J 

)lal0341-Il 

17.6 

16.2 

1.4 

10.0 

10.51 

1.54 

1.67 

2.65 

58.4 

41.6 

31.4 

13.2 

8.5 

7.8 

_J 

)lal3341-I2 

25.3 

22.6 

WSM 

13.0 

13.84 

1.63 

1.83 

2.65 

61.7 

38.3 

51.0 

19.5 

11.9 

10.7 

)lal  634 1-13 

29.1 

26.0 

3.1 

16.0 

17.17 

1.51 

1.69 

2.65 

57.2 

42.8 

42.3 

18.1 

11.9 

10.7 

)lal9341-I4 

36.0 

32.2 

3.8 

19.0 

20.50 

1.57 

1.76 

2.65 

59.3 

40.7 

46.0 

18.7 

11.9 

10.7 

_E 

Ial5511-J4 

30.0 

22.9 

7.1 

15.0 

15.77 

1.45 

2.05 

2.65 

54.9 

45.1 

99.9 

45.1 

31.1 

23.7 

_E 

Ial8511-J3 

37.2 

29.1 

8.1 

18.0 

19.07 

1.53 

2.05 

2.65 

57.7 

42.3 

99.9 

42.3 

27.7 

21.7 

_E 

Ia21511-J2 

44.7 

35.8 

8.8 

21.0 

22.37 

1.60 

2.05 

2.65 

60.5 

39.5 

100.1 

39.5 

24.7 

19.8 

J 

Ia24511-Jl 

53.1 

44.0 

9.0 

24.0 

25.67 

1.72 

2.05 

2.65 

64.9 

35.1 

100.0 

35.1 

20.5 

17.0 
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Appendix  Al.  Soil  conditions  at  time  of  acoustic  evaluation  (continued). 


Sample  File 
Number 


sacl6111-A4 


sacl9111-A3 


sac22111-A2 


sac25111-Al 


sacl6112-Bl 


sacl9112-B2 


sac22112-B3 


sac25112-B4 


sac08141-Cl 


saclll41-C2 


sacl4141-C3 


sacl6211-D4 


sacl9211-D3 


sac22211-D2 


sac25211-Dl 


sacl6212-El 


sacl9212-E2 


sac222l2-E3 


sac25212-E4 


sacl0241-Fl 


sacl3241-F2 


sacl6241-F3 


sacl9241-F4 


sacl33ll-G4 


sacl6311-G3 


sacl9311-G2 


sac22311-Gl 


sacl0341-Hl 


sacl3341-H2 


sacl6341-H3 


sacl7341-H4 


sacl4511-J7 


sacl6511-J3 


sacl7511-J6 


sacl9511-J2 


sac20511-J5 


sac23511-J4 


sac25511-Jl 


Tub 

Depth 

cm 


16.0  17.02 


19.0  20.35 


22.0  23.68 


25.0  27.01 


16.0  16.91 


19.0  20.21 


22.0  23.51 


25.0  26.81 


7.71 


11.0  11.59 


14.0  14.92 


17.0  18.25 


16.0  17.01 


19.0  20.34 


22.0  23.67 


25.0  27.00 


16.0  17.10 


19.0  20.43 


22.0 


25.0  127.09 


10.0  10.36 


13.0  13.69 


16.0  17.02 


19.0  20.35 


13.0  13.77 


16.0  17.10 


19.0  20.43 


22.0  23.76 


10.0  10.51 


13.0  13.84 


16.0  17.17 


17.0  18.28 


14.0  14.67 


16.0  16.87 


17.0  17.97 


19.0  20.17 


20.0  21.27 


23.0  24.57 


23.5  125.121 


Total 

Solid 

% 

Total 

Pores 

% 

HP 

Filled 

Pores 

% 

54.3 

45.7 

5.8 

51.6 

48.4 

5.2 

55.5 

44.5 

6.1 

57.2 

42.8 

6.5 

46.5 

53.5 

wm 

48.2 

51.8 

IKS 

48.3 

51.7 

Ea 

50.1 

49.9 

Kil 

50.5 

49.5 

5.0 

52.9 

47.1 

5.5 

53.9 

46.1 

mm 

54.8 

45.2 

mm 

38.5 

61.5 

18.6 

39.6 

60.4 

21.0 

40.3 

59.7 

21.6 

39.4 

60.6 

20.8 

33.9 

66.1 

15.2 

34.2 

65.8 

15.4 

35.1 

64.9 

16.0 

36.0 

64.0 

16.7 

47.9 

52.1 

28.1 

51.5 

48.5 

32.6 

51.3 

48.7 

32.2 

54.3 

45.7 

36.4 

46.4 

53.6 

50.5 

47.3 

52.7 

52.4 

46.2 

53.8 

50.6 

47.6 

52.4 

53.7 

63.0 

37.0 

88.9 

59.7 

40.3 

77.3 

60.2 

39.8 

78.8 

62.4 

37.6 

86.5 

50.9 

49.1 

100.1 

50.9 

49.1 

100.1 

50.9 

49.1 

100.1 

50.9 

49.1 

100.1 

50.9 

49.1 

100.1 

50.9 

49.1 

100.1 

50.9 

49.1 

100.1 
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Appendix  A2.  Soil  data  acquisition  system’s  autoexec.bat  file. 


@ECHO  OFF 

PROMPT  $_$t$_$d$_$p$g 

PATH 

C:\WINDOWS;C:\WINDOWS\COMMAND;C:\SyQuest;\WINDOWS;\DOS;Nxtgold;\ncsa2;\c60CKb 

in;Nc6(XNjinb\rich\Npkzip; 

PATH  C:\NU;%PATH% 

SET  SYMANTEC=C:\SYMANTEC 

SET  NU=C:\NU 

SETTEMP=C:\DOS 

SET  MOUSE=C:\MOUSE 

SET  INCLUDE=C:\C60(]MNCLUDE; 

SETLIB=C:\C600\LIB; 

SET  HELPFILES=C:\C600\HELP\*.HLP 
SET  INrr=C;\C600MNrr 
set  qhpib=c:Nhpib 
set  gpc=c:Ngraphic\fonts 

C;\WINDOW^COMMAND\MODE  CON:RATE=32  DELAY=1 

rem  -  By  Windows  Setup  -  LH  C:\WINDOWS\COMMAND\MSCDEX.EXE  /D:MTMIDE01 
/M:10 

rem  -  By  Windows  Setup  -  C;\MOUSE\MOUSE.EXE  /Q 

LH  C:\G095\DMS  MONITOR 

doskey 

c:\mouse\mouse.exe 

C:\WINDOWS^OMMAND\doskey 

ECHO. 

ECHO  TYPE  "WIN"  TO  START  WINDOWS  FOR  WORKGROUPS  3.1 1 

rem  cd  soil  " 

rem  C:\NU\NDD  C:/Q 

rem  C:\NUMMAGEC: 

c:Nncsa2Npkt8000  0x60  10  0x220  OxeOOO 

c:\winsock\winpkt  0x60 
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Appendix  A3.  Soil  data  acquisition  system’s  config.sys  file. 

rem  DEVICE=C:'^USLOGIC\BTDOSM.SYS  /D 

DEVICE=C:\DOS\HIMEM.SYS 

DEVICE=C:\LIB\SP0126.SYS 

DEVICE=C:^LUGPLAY\DRIVERS\DOS\DWCFGMG.SYS 

DEVICE=C:\DOSNEMM386.EXE  NOEMS  /X=C800-CBFF  /X=DC00-DFFF 

BIJFFERS=30 

FILES=45 

DOS=UMB 

DEVICEHIGH=C:MDOSNSETVER.EXE 

DOS=HIGH 

DEVICEfflGH=C:\MTM\MTMCDAI.SYS  /D:MTMIDE01  /P:  170,15 
SHELL=C:\DOS\COMMAND.COM  C;\DOS\/E:2048  /P 
DEVICE=C;\WINDOWS\IFSHLP.SYS 
STACKS=9,256 

rem  to  fix  error#12  to  the  EMM386  30.12.95 
remSTACKS=18,512 

remDEVICE=C:NSyQuestNSYQASPI\ASPI2DOS.SYS 

remDEVICE=C:\SyQuest\SCSI\SQDRIVER.SYS 

LASTDRr/E=Z 
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source  =  soil4 
cflags  =  /c  /AL  /Gs 
cdebug  =  /Zi  /Od 
Iflags  = 

Idebug  =  /co  /st:0x20(X) 

all:  $(source).exe  $(source).obj  msc_sub.obj 

$  (source). obj:  $(source).c  dtst_xmm.h  dtst_err.h  dtst_tls.h 
cl  $(cflags)  $(source).c 

msc_sub.obj:  msc_sub.c  dtst_xmm.h  dtst_err.h  dtst_tls.h 
cl  $(cflags)  msc_sub.c 

$(source).exe:  $(source).obj 

link  $(lflags)  $(source)„NUL,c:\c600MibMlibc7.1ib+clhpib  dtswtool 
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Appendix  A5.  Soil  data  acquisition  system’s  DTSTJTHS.H  file. 


Module:  DTST_TLS.H 
Function:  This  file  is  the  Data  Translation,  Inc. 
definition  file  for  the  software  tools 
called  from  Microsoft  C  V5. 1 


#defineAD_SECTION  0x0000 
#defineDA_SECTION  0x0001 
#define  AD_DA_SECTION  0x0002 


/*  Select  the  ad  section  */ 
/*  Select  the  da  section  */ 
/*  Select  both  sections 


SIMUL.  START  */ 


typedef  struct  BOARDS_STRUCT  { 

int  num_units;  /*  number  of  boards  installed  */ 

int  board_id;  /*  device  identifier  (see  manual)  */ 


base_address; 
driver_version  [10]; 
dma_channel_l ; 
dma_channel_2; 
interrupt_level; 
board_timeout; 
mux_id; 
ad_setup_bits; 
da_setup_bits; 
ad_sect_range; 
da_sect_range0; 
da_sect_rangel; 
chan_onbrd; 
chan_offbrd; 
ihrsh; 


}BOARDS; 


/*  base  address  of  this  board  */ 

/*  driver  name  asciiz  format  */ 

/*  1st  dma  channel  being  used  */ 

/*  2nd  dma  channel  being  used  */ 

/*  interrupt  being  used  */ 

/*  device  timeout  in  seconds  */ 

/*  multiplexer  id  (reserved)  */ 

/*  a/d  acquisition  setup  bits  */ 

/*  d/a  acquisition  setup  bits  */ 

/*  current  ad  input  range  in  millivolts  */ 

/*  current  da  0  input  range  in  millivolts  */ 

/*  current  da  1  input  range  in  millivolts  */ 

/*  on  board  input  channels  */ 

/*  off  board  input  channels  */ 

I*  dac  0  sample  for  theshold  triggering  */ 


/**  definitions  for  the  ad_setup_bits  field  of  the  BOARDS  structure  **/ 


#define  AD_DATA  0x4000 
#defineAD_POLARITY  0x2000 
#define  AD.INPUTS  Ox  1 000 
#define  AD_CLOCK  0x0800 
#defineAD_TSEL  0x0400 

#define  AD_TPOL  0x0200 

#defineAD_TRIGGER  0x0100 


/*  1  =  ad  binary  data,  0  =  ad  2s  comp  data  */ 
/*  1  =  ad  unipolar,  0  =  ad  bipolar  */ 

/*  1  =  se  inputs,  0  =  diff  inputs  */ 

/*  1  =  ext  clock,  0  =  int  clock  */ 

/*  1  =  analog  threshold  0  =  digital  */ 

/*  1  =  rising,  0  =  falling  */ 

/*  1  =  ext  clock,  0  =  int  clock  */ 


/**  definitions  for  the  da_setup_bits  field  of  the  BOARDS  structure  **/ 


#defmeDA_POLARl  0x0800 
#defineDA_POLAR0  0x0400 
#defineDA_DATAl  0x0200 
#defineDA_DATA0  0x0100 
#define  DA_TSEL  0x0080 

#defineDA_TPOL  0x0040 

#defmeDA_TRIGGER  0x0020 


/*  1  =  dac  1  unipolar,  0  =  da  bipolar  */ 

/*  1  =  dac  0  unipolar,  0  =  da  bipolar  */ 

/*  1  =  dac  1  binary  data,  0  =  da  2s  comp  data  */ 
/*  1  =  dac  0  binary  data,  0  =  da  2s  comp  data  */ 
/*  1  =  analog  threshold  0  =  digital  */ 

/*  1  =  rising,  0  =  falling  */ 

/*  1  =  ext  trigger,  0  =  int  trigger  */ 
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#defme  DA_CLOCK  0x0010  /*  1  =  ext  clock,  0  =  int  clock  */ 

#define  DAC_1  0x0002  /*  enable  dac  1  */ 

#define  DAC_0  0x0001  /*  enable  dac  0  */ 

#define  NO_DACS  0x0000  /*  disable  both  dacs  */ 

typ^ef  struct  DIGIT ALIO_STRUCT  { 

int  port;  /*  digital  io  port  to  operate  on  */ 

int  command;  /*  digital  io  command  */ 

int  direction;  /*  digital  io  direction  (input/output)  */ 

long  dio_value;  /*  digital  value  for  input  or  output  */ 

}  DIGITALIO; 

/**  definitions  for  the  port  field  of  the  DIGITALIO  structure  **/ 

#define  PORT_l  0x0002  /*  select  dio  port  1  */ 

#define  PORT_0  0x0001  /*  select  dio  port  0  */ 

/**  definitions  for  the  command  field  of  the  DIGITALIO  structure  **/ 

#define  SYNCH_MODE  0x0003  /*  select  synchronized  dio  mode  */ 

#define  SET_XFER  0x0002  /*  set  direction  and  then  set  /  get  value  */ 

#define  XFER_VAL  0x000 1  /*  set  /  get  a  port  value  */ 

#define  SET_DIR  0x0000  /*  set  port  direction  command  */ 

/**  definitions  for  the  direction  field  of  the  DIGITALIO  structure  **/ 

#define  DIO_OUTPUT  0x0001  /*  set  port  as  outputs  */ 

#define  DIO_INPUT  0x0000  /*  set  port  as  inputs  */ 

typedef  struct  FEATURES_STRUCT  { 

int  supported_features;  /*  major  subsystem  support  */ 

char  device_name[20];  /*  asciiz  encoded  device  name  */ 

char  module_name[20];  /*  asciiz  encoded  module  name  */ 

int  max_se_channels;  /*  max  number  of  ad  se  inputs  supported  */ 

int  max_di_channels;  /*  max  number  of  ad  di  inputs  supported  */ 

int  cg_list;  /*  depth  of  channel  gain  list  */ 

int  ad_support_bits;  /*  ad  support  flags  */ 

int  ad_sect_resolution;  /*  resolution  (in  bits)  of  ad  section  */ 

float  ad_sect_thruput;  /*  max  thruput  of  ad  section  */ 

float  ad_sect_clock;  /*  ad  section  base  clock  frequency  */ 

int  ad_sect_maxpos;  /*  max  positive  input  voltage  (in  millivolts)  */ 

int  ad_sect_maxneg;  /*  max  negative  input  voltage  (in  millivolts)  */ 

int  da_support_bits;  /*  da  support  flags  */ 

int  da_sect_number,  /*  number  of  dacs  on  board  */ 

int  da_sect_resolution;  /*  resolution  (in  bits)  of  da  section  */ 

float  da_sect_thruput;  /*  max  thruput  of  da  section  */ 

float  da_sect_clock;  /*  da  section  base  clock  frequency  */ 

int  da_sect_maxpos;  /*  max  positive  input  voltage  (in  millivolts)  */ 

int  da_sect_maxneg;  /*  max  negative  input  voltage  (in  millivolts)  */ 

int  digital_io_ports;  /*  number  of  digital  io  ports  */ 

int  digital_io_size;  /*  size  of  each  digital  io  port  */ 

int  num_user_ctrs;  /*  number  of  counter/timers  */ 

int  ex_support;  /*  extended  major  subsystem  support  flags  */ 
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}  BOARD_FEATURES; 

/**  definitions  for  the  supported_features  field  of  the  BOARD_FEATURES  structure  **/ 


#define  XSUP 
#define  SSH 
#define  TSCAN 
#define  AAF 
#define  ADFEFO 
#define  SCNFR 
#define  MUX 
#define  DTC 
#define  THTRIG 
#define  DSP 
#define  CAL 
#define  CLK 
#defineTMR 
#define  DIO 
#define  DA 
#define  AD 


0x8000 

0x4000 

0x2000 

0x1000 

0x0800 

0x0400 

0x0200 

0x0100 

0x0080 

0x0040 

0x0020 

0x0010 

0x0008 

0x0004 

0x0002 

0x0001 


/*  device  has  extended  support  */ 

/*  device  has  SSH  module  */ 

/*  device  supports  triggered  channel  scans  */ 

/*  device  has  anti-aliasing  filters  on  board  */ 

/*  device  has  FIFOs  on  input  data  stream  */ 

/*  device  supports  software  configurability  */ 

/*  device  supports  external  multiplexer  */ 

/*  device  supports  DT  Connect  */ 

/*  device  suppons  on  board  microprocessor  (intelligent) 
/*  device  supports  on  board  DSP  */ 

/*  device  supports  on  board  auto  calibration  */ 

/*  device  supports  programmable  clock  functions  */ 

/*  device  suppons  programmable  timer  functions  */ 

/*  device  supports  digital  input  /  output  */ 

/*  device  supports  da  */ 

/*  device  supports  ad  */ 


*/ 


!**  definitions  for  the  ad_support_bits  and  da_support_bits 

fields  of  the  BOARD_FEATURES  structure  **/ 


#define  PGL  0x0200 

#define  PGH  0x0100 

#define  SE  0x0080 

#define  DI  0x0040 

#define  BIN  0x0020 

#defineTWOS  0x0010 

#define  BIP  0x0008 

#define  UNI  0x0004 

#define  XCLK  0x0002 

#defineXTRIG  0x0001 


/*  ad  subsystem  has  PGL  module  */ 

/*  ad  subsystem  has  PGH  module  */ 

/*  a/d  supports  single  ended  inputs  */ 

/*  a/d  supports  differential  inputs  */ 

/*  adc  &  dac  support  binary  data  */ 

/*  adc  &  dac  support  twos  complement  data  */ 
/*  adc  &  dac  support  bipolar  inputs  */ 

/*  adc  &  dac  support  unipolar  inputs  */ 

/*  adc  &  dac  support  external  clock  */ 

/*  adc  &  dac  support  external  trigger  */ 


/**  definitions  for  the  ex_supp  field  of  the  BOARD_FEATURES  structure  **/ 


#defineTRIGLOUT  0x0004 
#define  SYNCHDO  0x0002 

#define  DT2896  0x0001 


/*  device  support  trigger  lockout  */ 

/*  device  supports  synch  dio  */ 

/*  device  supports  CGL  extension  thru  DT2896  */ 


typedef  struct  SET_ACQ_PARAMS_STRUC  { 


mt 

int 

float 

float 

float 


section; 

transfer_type; 

clock_rate; 

trig_rate; 

cut_off; 


)  SET_AO(i_PARAMS; 


/*  a/d  or  d/a  section  of  the  board  */ 
/*  type  of  acquisition  transfer  */ 
/*  acquisition  clock  rate  */ 

/*  re-trigger  rate  */ 

/*  cut  off  frequency  */ 


/**  definitions  for  the  transfer_type  field  of  the  SET_ACQ_P ARAMS  strucmre  **/ 

#define  BC_INTERRUPT  0x0010  /*  buffer  cycle  -  interrupt  driven  */ 

#define  BC_AUTOINIT  0x0011  /*  buffer  cycle  -  auto  init  dma  */ 
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#define  BJNTERRUPT  0x0020  /*  burst  -  interrupt  driven  */ 

#define  B_INTRPT_TRIG_SCAN  0x0024  /*  burst  -  interrupt  triggered  scan  */ 

#deflne  B_INTRPT_XTRIG_BUFF  0x0028  /*  burst  -  interrupt  xtrigger/buffer  */ 

#define  B_INTRPT_XTRIG_BUFF_TRIG_SCAN  0x002C  /*  burst  -  interrupt  xtriggered  scan 
*/ 

#define  B_SINGLE  0x0021  /*  burst  -  single  dma  channel  */ 

#define  B_SINGLE_TRIG_SCAN  0x0025  /*  burst  -  single  dma  triggered  scan  */ 

#define  B_SINGLE_XTRIG_BUFF  0x0029  /*  burst  -  single  dma  xtrigger/buffer  */ 

#define  B_SGL_XTRIG_BIJFF_TRIG_SCAN  0x002D  /*  burst  -  single  dma  xtriggered  scan 
*/ 

#define  B_DUAL  0x0022  /*  burst  -  dual  dma  channels  (gap  free)  */ 

#define  B_DUAL_TRIG_SCAN  0x0026  /*  burst  -  dual  dma  triggered  scan  */ 

#define  B_DUAL_XTRIG_BUFF  0x002A  /*  burst  -  dual  dma  xtrigger/buffer  */ 

#define  B_DUAL_XTRIG_BUFF_'1’RIG_SCAN  0x002E  /*  burst  -  dual  dma  xtriggered  scan  */ 

#define  C_INTERRUPT  0x0040  /*  continuous  -  interrupt  driven  */ 

#define  C_INTRPT_TRIG_SCAN  0x0044  /*  continuous  -  interrupt  triggered  scan  */ 

#define  C_INTRPT_XTRIG_BUFF  0x0048  /*  continuous  -  interrupt  xtrigger/buffer  */ 
#derine  C_INTRPT_XTRIG_BUFF_TRIG_SCAN  0x004C  /*  continuous  -  interrupt  xtriggered 
scan  */ 

#define  C_S INGLE  0x(X)41  /*  continous  -  single  dma  channel  */ 

#define  C_SINGLE_TRIG_SCAN  0x0045  /*  continuous  -  single  dma  triggered  scan  */ 

#define  C_SINGLE_XTRIG_BUFF  0x0049  /*  continuous  -  single  dma  xtrigger/buffer  */ 

#define  C_SINGLE_XTRIG_BUFF_TRIG_SCAN  0x004D  /*  continuous  -  single  dma 
xtriggered  scan  */ 

#define  C_DUAL  0x0042  /*  continous  -  dual  dma  channels  (gap  free)  */ 

#define  C_DUAL_TRIG_SCAN  0x0046  /*  continuous  -  dual  dma  triggered  scan  */ 

#define  C_DUAL_XTRIG_BUFF  0x004A  /*  continuous  -  dual  dma  xtrigger/buffer  */ 

#define  C_DUAL_XTRIG_BUFF_TRIG_SCAN  0x004E  /*  continuous  -  dual  dma  xtriggered 
scan  */ 

typedef  struct  SINGLE_STRUC  { 

int  sing_chan;  /*  channel  for  single  value  ad  operation  */ 

int  sing_gain;  /*  gain  for  single  value  ad  operation  */ 

int  sing_dacs;  /*  dac  mode  select  */ 

}  SINGLE; 

typedef  struct  TIME_DATE_STRUC  { 


int 

year; 

/*  year  of  time/date  stamp  */ 

int 

month; 

/*  month  of  time/date  stamp  */ 

int 

day; 

/*  day  of  time/date  stamp  */ 

int 

hour, 

/*  hour  of  time/date  stamp  */ 

int 

min; 

/*  minate  of  time/date  stamp  */ 

int 

sec; 

/*  seconds  of  time/date  stamp  */ 

}  T[ME_DATES; 

typedef  struct  COUNTER_TIMER_STRUC  { 

int 

ct_command; 

/*  counter/timer  command  */ 

int 

counter_select; 

/*  counter  to  operate  on  */ 

int 

clock_select; 

/*  0  =  internal;  1  =  external  */ 
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/*  countdown/divider  count  */ 
/*  users'  event  counter/flag  */ 
/*  measured  frequency  */ 


long  clock_count; 
long  event_flag; 
long  frequency; 

)  COUNTER_TIMERS; 

!***  ct_command  definitions  ***/ 

#define  RESET_COUNTER  0 
#define  SQUARE.WAVE  1 
#define  GENERA'IE_RATE  2 
#define  SW_ONE_SHOT  3 
#defineHW_ONE_SHOT  4 
#define  MEASURE_FREQ  5 
#define  COUNT_EVENT  6 
#derine  GET_COUNT  7 

/***  clock_select  definitions  ***/ 

#define  CT_INTERNAL  0 
#define  CT_EXTERNAL  1 

!***  counter_select  definitions  ***! 


/*  reset  the  selected  counter  */ 

/*  generate  a  square  wave  */ 

/*  generate  pulses  at  specified  rate  */ 
/*  software  triggered  one  shot  */ 

/*  hardware  triggered  one  shot  */ 

/*  measure  frequency  */ 

/*  count  events  *! 

I*  get  current  clock  count  */ 


/*  use  internal  clock  (600KHz)  */ 
/*  use  external  clock  */ 


#define  DT2812_TMR_CTR  0  /*  use  timer/counter  0  */ 


/***  channel  scan  enable  -  used  for  auto-channel  sequence  ***/ 

#define  CHANNEL_SCAN  0x80  /*  enable  channel  scan  */ 

/* . : . 

function  prototypes 

. */ 


int  far  cdecl  dt_counter_timer  (int  unit,  int  dev_handle, 

COUNTER_TIMERS  far  *cnttmrs); 
int  far  cdecl  dt_check_buffer  (int  unit,  int  section,  int  far  *buf_handle, 
int  far  *buf_status); 

int  far  cdecl  dt_create_buffer  (int  unit,  int  section,  int  far  *buf_size, 
int  far  *buf_ptr,  int  far  *buf_handle); 
int  far  cdecl  dt_create_cgl  (int  num_ch_ga,  int  far  *chan_array, 
int  far  *gain_array); 

int  far  cdecl  dt_delete_buffer  (int  unit,  int  section,  int  buf_handle); 
int  far  cdecl  dt_get_acq  (int  unit,  int  dev_handle,  int  section, 

SET_ACQ_PARAMS  far  *set_acq_params); 
int  far  cdecl  dt_get_board  (int  unit,  int  dev_handle, 

BOARDS  far  *brd_struct); 

int  far  cdecl  dt_get_buf_tstamp  (int  buffer_handle,  TIME_DATES  far  *tstamp); 
int  far  cdecl  dt.  get  features  (int  unit,  int  dev_handle, 

BOARD_FEATURES  far  *fea_struct); 
int  far  cdecl  dt_initialize  (char  far  *dev_name,  int  far  *dev_handle); 
int  far  cdecl  dt_link_xm_buffer  (int. unit,  int  section,  int  far  *xm_buffer_size, 
long  xm_buffer_ptr,  int  far  *buffer_handle); 
int  far  cdecl  dt_read_727_channel  (int  unit,  int  dev_handle,  int  channel, 
int  gain,  int  far  *data_ptr); 
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int  far  cdecl  dt_reset  (int  unit,  int  dev_handle); 

int  far  cdecl  dt_reset_buffer  (int  buffer_handle); 

int  far  cdecl  dt_set_acq  (int  unit,  int  dev_handle,  int  section, 

SET_ACQ_PARAMS  far  *set_acq_params); 
int  far  cdecl  dt_set_board  (int  unit,  int  dev_handle, 

BOARDS  far  *brd_struct); 
int  far  cdecl  dt_set_dio  (int  unit,  int  dev_handle, 

DIGITALIO  far  *dio_struct); 

int  far  cdecl  dt_single_acq  (int  unit,  int  dev_handle,  int  section, 

long  far  *sing_samples,  SINGLE  far  *sngl_operation); 
int  far  cdecl  dt_start_acq  (int  unit,  int  dev_handle,  int  section); 
int  far  cdecl  dt_stop_acq  (int  unit,  int  dev_handle,  int  section); 
int  far  cdecl  dt_tenninate  (int  dev_handle); 

int  far  cdecl  dt_wait_buffer  (int  unit,  int  section,  int  far  *buffer_handle); 
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Appendix  A6.  Soil  data  acquisition  system’s  DTST_XMM.H  file. 

y* _ * 

♦ _ * 

Module:  DTST_XMM.H  -* 

Function:  This  file  is  the  Data  Translation,  Inc.  -* 

definition  file  for  the  extended  memory  -* 

management  functions  using  Microsoft  -* 

HIMEM.SYS  called  from  Microsoft  C  V5.1.  -* 

_ * 

- */ 

int  far  cdecl  dt_xm_allocate_block  (int,  int  far  *); 
int  far  cdecl  dt_xm_free_block  (int) ; 

int  far  cdecl  dt_xm_get_emb_info  (int,  int  far  *,  int  far  *,  int  far  *); 

int  fai'  cdecl  dt_xm_get_free_mem  (int  far  *,  int  far  *); 

int  far  cdecl  dt_xm_^et_version  (int  far  *); 

int  far  cdecl  dt_xm_initialize  (); 

int  far  cdecl  dt_xm_lock_block  (int,  long  far  *); 

int  far  cdecl  dt_xm_move_dos_to_dos  (int  far  *,  int  far  *,  long); 

int  far  cdecl  dt_xm_move_dos_to_xm  (int  far  *,  int,  long,  long); 

int  far  cdecl  dt_xm_move_xm_to_dos  (int,  long,  int  far  *,  long); 

int  far  cdecl  dt_xm_move_xm_to_xm  (int,  long,  int,  long,  long); 

int  far  cdecl  dt_xm_reallocate_block  (int,  int); 

int  far  cdecl  dt_xm_unlock_block  (int); 


♦_ 

*. 

*. 

*. 

* 
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Appendix  A7.  Soil  data  acquisition  system’s  DTSTJERR.H  file. 


*  FILE  -  DTST_ERR.H  * 

*  DTT  Software  Tools  Error  Codes  * 

]|c  4;  4;  :ic  4c  4c  :{c  lie  :ic  4;  :ic  :<c  :ic  :ic  :t;  :jc  :|c  :|c  :|c  :ic  %  4: 4:  :)<  :)c %  4cl|<  *  %  4: 4c>|c  4c  1);  4c ;)( ;)c  3ic  :|c  4:  *  :j(  *  4: 4c  *  *  *  *  *  *  *  *  4c  iN  *  ^  ^  *  *  iic  :jc 


This  file  contains  the  DTI  Software  Tools  error  codes.  The  error 
codes  are  a  word  in  length.  The  upper  4  bits  encode  the  error  value 
into  one  of  16  different  range.  This  provides  for  a  count  of  4096 
different  error  codes  for  each  range. 


Range 

Type 

0x0000  -  OxOFFF 

INFORMATION  (not  necessarily  an  error) 

0x1000  -  OxlFFF 

AD  SECTION 

0x2000  -  0x2FFF 

DA  SECTION 

0x3000  -  0x3FFF 

BUFFERING 

0x4000  -  0x4FFF 

SYSTEM 

0x5000  -  0x5FFF 

DIGITAL  10  SECTION 

0x6000  -  0x6FFF 

EXTENDED  MEMORY 

0x7000  -  0x7FFF 

COUNTER/TIMER  SECTION 

0x8000  -  OxFFFF 

RESERVED 

/* . 

1 

INFORMATION  CODES 

4c/ 


#define  errComplete 
#defineerrInProcess 
#defineerrInvalidSection 
#defmeerrNoIOInProcess 
#define  errStoppedByUser 
#define  errIOInProcess 
#define  errDriverNotInit 


0x0000  /*  command  complete  */ 

0x0001  /*  command  in  process  */ 

0x0002  /*  invalid  board  section  specified  */ 

0x0003  /*  there  is  no  analog  i/o  in  process  */ 

0x0004  /*  user  stopped  function  in  process  */ 

0x0005  /*  analog  i/o  in  process  */ 

0x0006  /*  driver  has  not  been  initialized  */ 


#define  errToolkitAlreadyInit  0x0007  /*  toolkit  has  previously  been  initialized  */ 


/* 


AD  SECTION  CODES 


-*/ 


#define  errAdTrigger  Ox  1 000  /*  ad  trigger  error  */ 

#defineerrAdOvemin  0x1001  /*  ad  overrun  error  */ 

#defineerTAdNotReady  0x1002  /*  ad  system  not  ready  */ 

#defineeiTAdXferUnsupported  0x1003  /*  unsupported  ad  xfer  type  */ 

#defineerrInvalidChannel  0x1004  /*  invalid  channel  in  cgl  */ 
#defineerrInvalidGain  0x1005  /*  invalid  gain  in  cgl  */ 

#defineerrInvalidNumEntries  0x1006  /*  illegal  number  of  cgl  entries  */ 

#define  errlnvalidAdSetup  0x1007  /*  illegal  ad  setup  bits  in  SET_BOARD  */ 
#define  errlnvalid727 Channel  0x1008  /*  invalid  dt727  channel  selected  */ 

#define  errChannelSequence  0x1009  /*  illegal  channel  sequence  */ 
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/* 


DA  SECTION  CODES 


*/ 


I 


#defineeiTDaTrigger  0x2000  /*  da  trigger  error  */ 

^defineeirDaUnderrun  0x2001  /*  da  underrun  error  */ 

#defuieerrDaNotReady  0x2002  /*  da  system  not  ready  */ 

#defineerrDaXferUnsupported  0x2003  /*  unsupported  da  xfer  type  */ 

#defme  errIllegalDac  0x2004  /*  illegal  value  for  dac  selected  */ 

#define  errInvalidDaSetup  0x2005  /*  illegal  da  setup  bits  in  SET_BOARD  */ 


_ _ _ _ _ _ _ 

I  BUFFERING  CODES 

-  */ 


#defineeiTBufferNotFound  0x3000  /*  Buffer  is  not  in  the  BTL  */ 

#defmeerrBufferSize  0x3001  /*  illegd  buffer  size  requested  */ 

#defineerrNoMoreBCBs  0x3002  /*  there  are  no  more  BCB's  available  */ 

#define  errBufferOverrun  0x3003  /*  buffer  not  returned  in  time  */ 

#define  errBufferNotReady  0x3004  /*  buffer  not  ready  for  io  */ 

#define  errBuffOnOddAddress  0x3005  /*  buffers  must  reside  on  a  word  boundry  */ 
#define  errBuffXsDmaBoundry  0x3006  /*  Buffer  croses  dma  page  boundry  */ 
#define  errNoBuffers  0x3007  /*  There  are  not  any  BCB's  in  the  BTL  */ 

#define  errElegalBtl  0x3008  /*  Illegal  BTL  for  requested  operation  */ 


/* 


SYSTEM  CODES 


■*! 


#defmeerrUnknownCommand  0x4000  /*  unknown  command  error  */ 

#defineerrInvalidLun  0x4001  /*  invalid  logical  unit  error  */ 

#defineerrTimeout  0x4002  /*  timeout  error  *! 

#defineerrDmaConflict  0x4003  /*  dma  channel  conflict  error  */ 

#defmeerrInterruptConflict  0x4004  /*  interrupt  channel  conflict  error  */ 

#defineerrIllegaIDma  0x4005  /*  illegal  dma  channel  selection  */ 

#defineerrlllegallnterrupt  0x4006  /*  illegal  interrupt  channel  selection  */ 
#defineerrUnsupportedFunction  0x4007  /*  unsupponed  function  request  */ 

#defmeerrIllegalDriverCall  0x4008  /*  illegd  driver  level  call  detected  *! 
#defmeerrniegalBufferType  0x4009  /*  illegal  buffer  for  requested  op.  */ 

#defmeerrNoInterruptChannel  0x400A  /*  function  requires  int  chan  selection  */ 

#defineeiTNoDmaChannel  0x400B  /♦  function  requires  dma  chan  selection  */ 

#define  errlllegalRate  0x400C  /*  invalid  clock  rate  entered  */ 

#define  errIllegalScanRate  0x400D  /*  illegal  channel  scan  rate  */ 

#defme  errInvalidScanRate  0x400E  /*  invalid  scan  rate  for  selected  conversion  rate  */ 


^define  errReducedRate 
#define  errDrvNotFound 
#define  errTooMany  Files 
#defineerrInvalidAccess 
#define  errOpenOther 
#defme  errin  validHandle 
#defineerrCloseC)ther 
#defineerrInvalidFunction 
#defme  errSendOther 
#define  err727notAttached 


0x400F  /*  gain  selection  has  reduced  conversion  rate  */ 

0x4010  J*  driver  was  not  found  on  open  */ 

0x401 1  /*  cannot  open  driver,  too  many  files  open  */ 

0x4012  /*  cannot  open  driver,  invalid  access  specified  */ 
0x4013  /*  unknown  error  on  open  driver  call  */ 

0x4014  /*  invalid  device  handle  */ 

0x4015  /*  unknown  error  closing  driver  */ 

0x4016  /*  invalid  driver  function  requested  */ 

0x4017  /*  unknown  error  sending  driver  packet  */ 

0x4018  /*  no  dt727  connected  to  select^  unit  */ 
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#defme  errIllegalCutOff  0x4019  /*  illegal  input  cut-off  frequency  */ 

#define  errMinimumRate  0x40 1 A  /*  minimum  conversion  rate  not  met  */ 

#defme  errInvalidTimeout  0x40 1 B  /*  invalid  timeout  value  */ 

#define  errInvalidBtl  0x40 1C  /*  invalid  BTL  */ 

#define  errInvalidCgl  0x40 ID  /*  invalid  CGL  */ 


I  DIGITAL  10  SECTION  CODES  I 

- *! 

#defmeerrIllegalDioPort  0x5000  /*  illegal  digital  port  selection  */ 
#defmeerrIllegalDioCommand  0x5001  /*  illegal  distal  port  command  */ 

#defmeerrDioDirection  0x5002  /*  inv^id  digital  i/o  direction  */ 


/* 

i 


EXTENDED  MEMORY 


#defineerrXMSDrvrNotInstalled  0x6000  /*  XMS  driver  is  not  installed  */ 

#defineerTFuncNotImp  0x6001  /*  XMS  function  is  not  implemented  */ 

#defmeerrVDiskFound  0x6002  /*  VDisk  device  installed  */ 

#defineerrAllXMAllocated  0x6003  /*  all  extended  memory  is  allocated  */ 

#defineerrAllXMHandlesAlloc  0x6004  /*  all  XMS  handles  have  been  allocated  */ 

#defineerrInvalidXMSHandle  0x6005  /*  invalid  XMS  handle  used  */ 

#defineerrXMSBlockLocked  0x6006  /*  block  must  be  unlocked  for  selected  operation  */ 

#defineeiTOverflowLockCount  0x6007  /*  the  block  lock  count  has  overflowed  */ 

#defineeiTLockFailed  0x6008  /*  block  could  not  be  locked  */ 

#deftneerrBlockUnlocked  0x6009  /*  block  is  already  unlocked  */ 

#derineerrA20LineError  0x600A  /*  an  error  on  the  A20  address  line  has  occured  */ 

#defineerrInvalidSrcHdl  0x600B  /*  invalid  source  handle  */ 

#defmeerrInvalidSrcOff  0x600C  /*  invalid  source  offset 

#defineeiTlnvalidDestHdl  0x600D  /*  invalid  destination  handle  */ 

#defmeerrInvalidDestOff  0x600E  /*  invalid  destination  handle  */ 

#defmeerrInvalidLength  0x600F  /*  invalid  transfer  count  in  move  block  */ 
#defmeerrInvalidOverlap  0x6010  I*  cannot  perform  move  block,  invalid  block  overlap  */ 
#defineerrParity  0x6011  /*  cannot  perform  move  block,  parity  error  */ 

#defme  errlUegalBufSize  0x6012  /*  requested  buffer  size  is  invalid  */ 
#defineerrUnlmownXMSError  0x6013  /*  an  unknown  XMS  error  has  occured  */ 


/* 


COUNTERyriMER  SECTION 

*! 


#define  errInvalidCTCommand  0x7000  /*  Invalid  counter/timer  command  */ 
^define  errInvalidCTSelect  0x7001  /*  illegal  counter/timer  selected  */ 
#define  eirlnvalidClkSelect  0x7002  /*  illegal  clock  selected  */ 

^define  errInvalidClkCount  0x7003  /*  illegal  clock  count  selected  */ 

#defme  errInvalidSource  0x7004  /*  invalid  count  source  specified  */ 
#deflne  errInvalidGate  0x7005  /*  invalid  gate  specified  */ 

#define  errCountOverFlow  0x7006  /*  counter  has  overflowed  */ 

#define  errlnvalidPeriod  0x7007  /*  invalid  period  specified  */ 

#define  errInvalidRatio  0x7008  /*  invalid  ratio  specified  */ 
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Appendix  A8.  Soil  data  acquisition  system’s  MSC_SUB.C  file. 

^♦♦*****>|t*iNiic=l<=it******>l=***********************************************°*'*** 

Copyright  (C)  1990.  Data  Translation,  Inc.,  100  Locke  Drive, 

Marlboro  Massachusetts  01752-1 192. 

All  rights  reserved.  This  software  is  furnished  to  purchaser  under  a 
license  for  use  on  a  single  computer  system  and  can  be  copied  (with 
the  inclusion  of  DTI's  copyright  notice)  only  for  use  in  such  system, 
except  as  may  be  otherwise  provided  in  writing  by  Data  Translation, 

Inc.,  100  Locke  Drive,  Marlboro,  Massachusetts  01752-1 192. 

The  information  in  this  document  is  subject  to  change  without  notice 
and  should  not  be  construed  as  a  commitment  by  Data  Translation,  Inc. 

Data  Translation,  Inc.  assumes  no  responsibility  for  any  errors  that 
may  appear  in  this  document. 

Data  Translation  cannot  assume  any  responsibility  for  the  use  of  any 
portion  of  this  software  on  any  equipment  not  supplied  by  Data 
Translation,  Inc. 

j(e  *  :^e  *  3^  :ie  *  :ic  *  lie  ^  if:  :|c  *  :if<  *****  3}«  *♦*  ^  **’>'********************  */ 

#include  <stdio.h> 

#include  <conio.h> 

#include  "dtst_tls.h" 

#include  "dtst_err.h" 

#include  "dtst_xmm.h" 

void  wait_for_user  0 

{ 

printf("\ri\n">  Keystroke  to  continue..."); 
flushallO; 
getcharO; 

} 

void  featur  ( int  unit,  int  dev_handle) 

{ 

BOARD_FEATURES  features; 
int  status; 

printf  ("\nGetting  Device  Features:  STATUS  =  0x%x\n", 

status  =  dt_get_features  (unit,  dev_handle,  &features)); 

if  (status  ==  errComplete) 

printf  ("NnFeatures  Word . Ox%x",features.supported_features); 

printf  ("\nDevice  Name . %s"  ,features.device_name); 

printf  ("\nModule  Name . %s"  ,features.module_name); 

printf  ("\nlnput  SE  Channels . %d''  ,features.max_se_channels); 

printf  ("\nlnput  DI  Channels . %d"  ,features.max_di_channels); 

printf  ("\nCGL  Size . %d"  ,features.cg_list); 

printf  ("\nAD  Section  Suppon . Ox%x'',features.ad_support_bits); 

printf  ("\nAD  Resolution . %d"  ,features.ad_sect_resolution); 

printf  ('VAD  Thruput . %.2f',features.ad_sect_thruput); 
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printf  ("\nAD  Base  Clock . %.2f',features.ad_sect_clock); 

printf  ("\nAD  Max  +  Input . %d"  .features. ad_sect_maxpos); 

printf  ("\nAD  Max  -  Input . %d"  .features. ad_sect_maxneg); 

printf  ("\nDA  Section  Support . Ox%x".features.da_support_bits); 

printf  ("\nNumber  DACS . %d"  .features.da_sect_number); 

printf  ("\nDA  Resolution . %d"  .features.da_sect_resolution); 

printf  ("\nDA  Thruput . %.2f'.features.da_sect_thruput); 

printf  ("NnDA  Base  Clock . %.2f'.features.da_sect_clock); 

printf  ("\nDA  Max  +  Output . %d"  .features.da_sect_maxpos); 

printf  ("XnDA  Max  -  Output . %d"  .features.da_sect_maxneg); 

printf  ("ViNumber  DIO  Ports . %d"  .features.digital_io_ports); 

printf  ("\nSize  DIO  Ports  (bits)....%d"  .features.digital_io_size); 
printf  ("SnNum  user  counter/timers..%d"  .features.num_user_ctrs); 
printf  ('^Extended  support . %d"  .features.ex_support); 


wait_for_user(); 

}  /*  end  featur  routine  */ 


%  4c  life  :fc  % 4c  3k  %  *  3f:  %  %  :f:  sK  %  ^  ^  ^  :4c  :fc  :fc  %  %  sfc  :if:  }}c  %  %  3k  3k ’fe  ^  ^ 

function  put_boar 


This  routine  tests  the  dt_set_board  toolkit  routine. 

void  put_boar  (int  unit,  int  dev_handle,  int  dmal,  int  dma2, 
int  intrpt,  int  ad_setup,  int  da_setup,  int  tmout) 

{ 

BOARDS  board; 

board.dma_channel_l  =dmal; 
board.dma_channel_2  =  dma2; 
board.interruptjevel  =  intrpt; 
board.board_timeout  =  tmout; 
board.ad_setup_bits  =  ad_setup; 
board.da_setup_bits  =  da_setup; 
board.chan_onbrd  =  0; 
board.chan_offbrd  =  0; 
board.thrsh  =  0; 


printf  ("\nSetting  device's  BOARD  structure:  STATUS  =  0x%x\n", 
dt_set_board  (unit,  dev_handle,  &board)); 

}  /*  end  routine  put_boar  */ 


function  get_boar 
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This  routine  tests  the  dt_get_board  toolkit  routine. 

♦  ♦********)|c**!(c*****:(c*****s|c  ********  Jlcic^********************************^ 

void  get_boar  ( int  unit,  int  dev_handle ) 


BOARDS  board; 
int  status; 

printf  ("\nGetting  device's  BOARD  structure  STATUS:  =  Ox%x\n", 
status  =  dt_get_board  (unit,  dev_handle,  &board)); 

if  (status  ==  eirComplete) 

{ 

printf  ("\nNumber  Installed  Units...%d"  ,board.num_units); 

printf  ("\nDevice  ID . Ox%x"  ,board.board_id); 

printf  ("SnDevice  Base  Address . Ox%x"  .board. base_address); 

printf  ('ViDevice  Driver  Version... .%s"  , board. driver_version); 

printf  ("\nDMA  Channel  1 . %d"  ,board.dma_channel_l); 

printf  ("\nDMA  Channel  2 . %d"  ,board.dma_channel_2); 

printf  ("\nlnterrupt  Channel . %d"  ,board.interrupt_level); 

printf  ("ViTimeout  Value  (secs) . %d"  ,board.board_timeout); 

printf  ("\nA/D  Setup  Value . 0x%04x", board. ad_setup_bits); 

printf  ("\nD/A  Setup  Value . 0x%04x", board. da_setup_bhs); 

printf  ('ViAD  Range . %d"  ,board.ad_sect_range); 

printf  ("\nDA0  Range . %d"  , board. da_sect_rangeO); 

printf  ("\nDAl  Range . %d"  ,board.da_sect_rangel); 

printf  ("\nMuxID . %2X"  ,board.mux_id); 

printf  ("\nChannels  onboard . %d"  ,board.chan_onbrd); 

printf  ('NnChannels  offboard . %d"  ,board.chan_offbrd);  ' 

printf  ("\nThreshold  trigger  value.. %d’'  .board. thrsh); 

} 

wait_for_user(); 

}  /*  end  routine  get_boar  */ 


^*  *************************************** 

function  get_acq 


This  routine  tests  the  dt_get_acq  toolkit  routine. 

void  get_acq  (int  unit,  int  dev_handle,  int  section) 


SET_Aai.PARAMS 
int  status; 


printf  ("VnGetting  acquisition  information  STATUS:  Ox%xVi", 
status  =  dt_get_acq  (unit,  dev_handle,  section,  &sap)); 


Apjjendices;  Page  A23 


Final  Report 


if  (status  ==  errComplete) 

{ 

printf  ("SnSection . 

printf  ("NnTransfer  Type . 

printf  ("\nConversion  Rate.. 

printf  ("ViTrigger  Rate . 

printf  ('NnCut  off  frequency, 

1 


%d"  ,sap. section); 

..Ox%x",sap.transfer_type); 

....%.2f',sap.clock_rate); 

.%.2f',sap.trig_rate); 

....%.2f',sap.cut_off); 


wait_for_user(); 


}  /*  end  get_acq  routine  */ 


**♦***********************************************************+***** 


function  set_acq 

This  routine  tests  the  dt_set_acq  toolkit  routine. 

:4e  rfc  ^  %  lie  4:  4;  3}; :)( ife  %  ^  )|c  :f;  %  %  ^  ^  9k  *  *  %  %  sfc  :jc  3)e  jfe  ^  9f:  ^  :4c  :4c  :f:  4: 


void  set_acq  (int  unit,  int  dev_handle,  int  section, 

int  tran_type,  float  rate,  float  trig_rate,  float  cut_off) 


{ 

SET_ACX^PARAMS  sap; 

sap.section  =  section; 
sap.transfer_type  =  tran_type; 
sap.clock_rate  =  rate; 
sap.trig_rate  =  trig_rate; 
sap.cut_off  =  cut_off; 

printf  ("\nSetting  Acquisition  Information  STATUS  =  Ox%x\n", 
dt_set_acq  (unit,  dev_handle,  section,  &sap)); 

}  /*  end  set_acq  routine  */ 
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Appendix  A9.  Soil  data  acquisition  system’s  soil4.c  file. 


soil4.c 

I  Written  for  U.S.  Army  CERL  Soil  Project 

I  Version  1.0  16.11.95 

I  1.1  12.12.95  Write  lots  o' data. 

I  1.2  14.12.95  Checks  for  bad  data  and  waits  for  I 

good. 

I  1.3  18.12.95  Includes  writes  to  comment  files. 

I  Close  to  user  version. 

I  1.4  27.12.95  Save  one  or  multiple  files. 

I  Ready  for  use. 

I  Written  by  Rich  Czerwinski  &  Nadine  Barrie  Smith 

I  &  Copyright  (C)  1990.  Data  Translation,  Inc.  & 

I  Functions: 

1  change_settings(char  *codes) 

I  collect(numchan,  numbuf,  numtrials,  bsize,  clk_rate) 

1  control_panel(int  *rpt,  int  *bur,  int  int  *amp) 

!  enter_data_file() 

I  enter_one_file() 

I  error_handler  (error,  routine) 

I  refresh_traces(int  rpt) 

I  sample(int  numchan,  int  numbuf,  int  bsize,  int  clk_rate) 

I  timerO 

I 

#include  <cfunc.h> 

#include  <chpib.h> 

#include  <sicl.h> 

#include  <string.h> 

#include  <stdio.h> 

#include  <stdlib.h> 

#include  <graph.h> 

#include  "dtst_tls.h" 

#include  ''dtst_err.h" 

#include  "dtst_xmm.h" 

#include  <sys\types.h> 

#include  <sysNtimeb.h> 

#define  NUMCHAN  4 
#define  NUMBUF  2 

#define  SAMPLES  2048 

#defineCLK_RATE  100000. 

#define  SANITYNUMTRIALS  1 
#define  NUMTRIALS  1 
#defineMAXAMP  5. 

#defineFl  59 
#defineF2  60 
#defineF3  61 
#defineF4  62 
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#defineF5  63 
#defineF6  64 
#defineF7  65 
#define  F8  66 
#define  ESC  27 

/* - 

I  Global  File  definations  I 

- */ 

FILE  *out; 

int 

j. 

real_trig, 

sampleO, 

change_settings(), 

buffer[NUMBUF]  [SAMPLES  *NUMCHAN] , 
isc  =  7, 

wait_time  =  2,  /*  wait  2s  before  next  action  */ 

source  =  716;  /*  HPIB  parameters 

*1 

void 

collectO, 

controLpanelQ, 

reffesh_tracesO; 

float 

localbuf[S  AMPLES  *NUMCHAN], 
data[NUMCHAN]  [SAMPLES] ; 

struct  timeb  xtime;  /*  for  the  timer  function  */ 

long 

timel, 

time2; 

double 

limits[4], 

minuslimits[4]; 

char 

str_buffer[100], 

file_name[20], 

name[50]; 

FILE 

*fp, 

*fopen(); 

mainO 

{ 


char 
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♦codes; 


int 

error, 
rpt  =  70, 
bur  =  5, 
frq  =  2000, 
amp  =  100, 

i, 

window_length  =  SAMPLES; 

limits[0]  =  5.; 
limits[l]  =  5.; 
limits[2]  =  3.5; 
limits[3]  =  1.25; 
mjiiuslimits[0]  =  -5; 
minuslimitsfl]  =  -5; 
minuslimits[2]  =-1.5; 
minuslimits[3]  =  0.25; 

reaLtrig  =  (int)(0.5  +0.5*rpt*CLK_RATE/(1000*NUMCHAN)); 


/* . 

I  initialize  HPIB 

. */ 

error  =  lORESET  (isc); 
error_handler  (error,  "lORESET"); 
error  =  lOTIlVEOUT  (isc,  5.0); 
eiTor_handler  (error,  "lOTlMEOUT"); 
error  =  lOCLEAR  (isc); 
error_handler  (error,  "lOCLEAR”); 

codes  =  "M7,CT0,T0,W1,H0,A0,C0,DTY  50  %,OFS  0.0, LO, FRQ  2  KHZ,  AMP 
100MV,BUR  5  #,  RPT  70  MS,  DO"; 
error  =  change_settings(codes); 
eiTor_handler  (error,  "lOOUTPUTS"); 

/* 

1  set  up  graphics  for  wave  display 

. .  */ 

_setvideomode(  _MAXRESMODE ); 

_clearscreen(  _GCLEARSCREEN ); 

for  (j  =  0;  j  <  NUMCHAN;  j-t-+) 

{ 

_setcolor(9); 

_setviewport(  1 0, 10+j  *70,600,75+j*70 ); 

_setwindow(0,  real_trig,  limits[j],  windowjength  -real_trig,minuslimits[j]); 
_rectangle_w(lGBORDER,real_trig,limitsIj],window_length- 
real_trig,minusliniitsU]); 

} 

control_panel(&rpt,  &bur,  &frq,  &amp); 

/*  if  ESC  is  selected  */ 
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strcpy(codes,  "FRQ  2000  HZ,  AMP  10  MV"); 
error  =  change_settings(codes); 
error_handler  (error,  "lOOUTPUTS"); 
_setvideomode(  _DEFAULTMODE ); 

}  /*  END  MAIN*/ 


I  Function:  refresh_traces(int  rpt) 

1 

- */ 

void  refresh_traces(int  rpt) 

{ 

int  i,  j,  windowjength  =  SAMPLES,  real_trig; 

real.trig  =  (int)(0.5  +  0.5*rpt*CLK_RATE/(1000*NUMCHAN)); 

!*- 

I  Clear  out  buffers 

. *! 

for(i  =  0;  i  <  NUMBUF;  i++) 

for  (j  =  0;  j  <  NUMCHAN*SAMPLES;  j  ++) 
buffer[i][j]  =  0; 

for  (j  =  0;  j  <  NUMCHAN*SAMPLES;  j  ++)  localbuflj]  =  0.; 

/* . 

I  set  up  graphics  for  wave  display 

. *! 

_setvideomode(  _MAXRESMODE ); 

_clearscreen(_GCLEARSCREEN ); 

/* 

I  Get  data  from  board 
- */ 

collect(NUMCHAN,  NUMBUF,  SANTTYNUMTRIALS,  NUMCHAN*SAMPLES, 
CLK.RATC); 

for  (j  =  0;  j  <  NUMCHAN;  j++) 

{ 

_setcoIor(9); 

_setviewport(10,10+j*70,600,75+j*70 ); 

_setwindow(0,  real_trig,  limits[j],  windowjength  -real_trig,minuslimits[j’]); 
_iectangle_w(_GBORDER,real_trig,limitsIj],window_length- 
real_trig,minuslimits[j] ); 

_setcolor(0); 

_moveto_w(  0,  -MAXAMP*data[j][real_trig]); 
for  (i  =real_trig;  i<SAMPLES;  i++) 

{ 

data[j][i]=MAXAMP*localbufli*NUMCHAN+j]; 

_setcolor(10+j); 
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)/* 


_lineto_w(  i,  -data[j][i] ); 

} 


} 

end  refTesh_traces  */ 


/* - 

I  Function:change_settings(char  *ccxies) 

I 

- */ 

int  change_settings(char  *codes) 

{ 

return  (lOOUTPUTS  (source,  codes,  strlen(codes) )); 

}  /*  end  change_settings  */ 


I* 

I 

I  Function:  controLpanel 
I 

. . . *! 

void  control_panel(int  *rpt,  int  *bur,  int  *frq,  int  *amp) 

{ 

static  char  codes[80]=  *rptstr,  *burstr,  *frqstr,  *ampstr,  key; 

static  int  error,  i,  enable  =  1; 

int  stop.freq; 

int  freq=0; 

int  amplitudes  10; 

int  x,y; 

char  freq_ext[5]; 
char  amp_ext[5]; 

while  (key  !=  ESC ) 

{ 

_settextcolor(9); 

_settextposition(  20, 1 ); 
sprintf(str_buffer,"[Fl]  FRQ  %d  Hz  ",  *frq); 

_outtext(  str_buffer ); 

_settextposition(  21, 1 ); 

sprintf(str_buffer,"[F2]  BUR  %d  cycles  per  burst  ",  *bur); 
_outtext(  str_buffer ); 

_settextposition(  22,  1 ); 
sprintf(str_buffer,"[F3]  RPT  %d  ms  ",  *rpt); 

_outtext(  str_buffer ); 

_settextposition(  23, 1 ); 

sprintf(str_buffer,"[F4]  AMP  %d  mV  ",  *amp); 

_outtext(  str_buffer ); 

_settextposition(  24, 1 ); 
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if  (enable) 

sprintf(str_buffer,"[F5]  Turn  that  !&*$@  thing  off!"); 

else 

sprintf(str_buffer,"[F5]  OK,  turn  it  back  on  again."); 
_outtext(  str_buffer ); 

_settextposition(  25, 1 ); 
sprintf(str_buffer."[F6]  Refresh  trace  display"); 

_outtext(  str_buffer ); 

_settextposition(  26, 1 ); 
sprintf(str_buffer,"lT7]  Save  to  Multiple  files"); 

_outtext(  str_buffer ); 

_settextposition(  20, 40 ); 
sprintf(str_buffer,"[F8]  Save  current  data  to  file"); 

_outtext(  str_buffer ); 

_settextposition(  27,  1  ); 
sprintf(str_buffer,"[ESC]  Quit"); 

_outtext(  str_buffer ); 


key  =  getchO; 
switch  (key) 

{ 

case  FI: 

_settextposition(  27, 1 ); 
sprintf(str_^buffer,"New  Frequency  (in  Hz):"); 
_outtext(  str_buffer ); 
scanf("%s",  frqstr); 

*frq  =  atoi(frqstr); 

strcat(codes,  "FRQ "); 

strcat(codes,  frqstr); 

strcat(codes,  "HZ"); 

error  =  change_settings(codes); 

eiTor_handler  (error,  "lOOUTPUTS"); 

_settextposition(  27, 1 ); 

sprintf(str_buffer,"  "); 

_outtext(  str_buffer ); 

control_panel(rpt,  bur,  frq,  amp); 

break;  /*  break  case  FI  */ 


case  F2; 

_settextposition(  27, 1 ); 

sprintf(str_buffer,"New  burst  length  (in  cycles):"); 
_outtext(  str_buffer ); 
scanf("%s",  burstr); 

*bur  =  atoi  (burstr); 

strcat(codes,  "BUR  "); 

strcat(codes,  burstr); 

strcat(codes, "  #"); 

error  =  change_settings(codes); 

error.handler  (eiror,  "lOOUTPUTS"); 
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_settextposiiion(  27, 1 ); 
sprintf(str_buffer,"  "); 

_outtext(  str_buffer ); 
control_panel(rpt,  bur,  frq,  amp); 
break;/*  break  case  F2  */ 


case  F3: 

_settextposition(  27,  1 ); 

sprintf(str_b'jffer,*'New  pulse  repitition  interval  (in  ms):"); 
_outtext(  str_buffer ); 
scanf("%s",  rptstr); 

*rpt  =  atoi(rptstr); 

strcat(codes,  "RPT  "); 

strcat(codes,  rptstr); 

strcat(codes,  "MS"); 

error  =  change_settings(codes); 

error_handler  (error,  "lOOUTPUTS"); 

_settextposition(  27,  1 ); 

sprintf(str_buffer,"  "); 

_outtext(  str_buffer ); 
control_panel(rpt,  bur,  frq,  amp); 
break;  /*  break  case  F3  */ 

case  F4: 

_settextposition(  27,  1  ); 

sprintf(str_buffer,"New  pulse  amplitude  (in  mV):"); 
_outtext(  str_buffer ); 
scanf("%s",  ampstr); 

*amp  =  atoi(ampstr); 
strcpy(codes,  "AMP "); 
strcat(codes,  ampstr); 
strcat(codes,  "MV"); 
error  =  change_settings(codes); 
eiTor_handler  (error,  "lOOUTPUTS"); 

_settextposition(  27,  1  ); 
sprintf(str_buffer,"  "); 

_outtext(  str_buffer ); 
control_panel(rpt,  bur,  frq,  amp); 
break;  /*  break  case  F4  */ 

case  F5: 

if  (enable) 

{ 

strcpy(codes,  "AMP  10  MV"); 
enable  =  0; 

] 

else 

{ 

strcpy(codes,  "AMP  "); 
itoa(*amp,  ampstr,  10); 
strcat(codes,  ampstr); 
strcat(codes, "  MV"); 
enable  =  1; 

) 
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error  =  change_settings(codes); 
error_haridler  (error,  "lOOUTPUTS"); 
control_paneI(rpt,  bur,  frq,  amp); 
break;  /*  break  case  F5  */ 


case  F6: 

refresh_traces(*rpt); 
control_panel(rpt,  bur,  frq,  amp); 
break; 


caseF7: 

en  ter_data_fileO; 
amplitude  =  5; 

stop_freq=37;  /*  stop  at  10  kHz  */ 

for  ( y  =  0;  y  <  5;  y++) 

{ 

amplitude  =  2  *  amplitude; 
itoa(y,  amp_ext,  10); 
fTeq=750; 

if  (  y  ==  3  II  y  ==  4  )  /*  stop_freq=6kHz  */ 

{ 

stop_freq  =  21; 

} 

for  (x  =  0;  X  <  stop_,freq;  x++) 

{ 

itoa(amplitude,  ampstr,  10); 
strcpy(codes,  "AMP "); 
strcat(codes,  ampstr); 
strcat(codes, "  MV, "); 

itoa(x,  freq_ext,  10); 

freq=ffeq+250; 
itoa(freq,  frqstr,  10); 
strcat(codes,  "FRQ "); 
strcat(codes,  frqstr); 
strcat(codes,  "  HZ"); 

_settextposition(  22,  5); 
sprintf(str_buffer,"  %s  ",  codes  ); 
_outtext(  str_buffer ); 

error  =  change_settings(codes); 
error.handler  (error,  "lOOUTPUTS"); 

strcpy(name,  "C:\NsoiI\\lata\\"); 
strcat(name,  file_name); 
strcat(name, "."); 
strcat(name,  amp_ext); 
if  (x  <  10)  strcat  (name,  "0"); 
strcat(name,  freq_ext); 

_settextposition(  23,  5); 
sprintf(str_buffer,"  %s  ",  name  ); 
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.outtext(  str_buffer ); 


timerO; 

fp  =  fopen(name,  "w"); 
collect(NUMCHAN,  NUMBUF, 
NUMTRIALS,  NUMCHAN*SAMPLES,  CLK_RATE); 

for  (i  =  real_trig;  i<SAMPLES;  i-H-) 

{ 

/*  fprintf(fp,  "%d  */ 

for  (j  =  0;  j  <  NUMCHAN;  j++) 

{ 

fprintf(fp,  "%f  ", 

MAXAMP*localbufri*NUMCHAN  +j]); 

) 

fprintf(fp,  "\n"); 

) 

fclose(fp); 

refresh_traces(*rpt); 

}  /*  end  for  (x  =  0;  x  <  20;  x++)  */ 

}  /*  end  for  (y  =  0;  x  <  6;  x++)  */ 

strcpy (codes,  "FRQ  2000  HZ,  AMP  10  MV"); 
error  =  change_settings(codes); 
erTor_handler  (error,  "lOOUTPUTS"); 

_setvideomode(  _DEFAULTMODE ); 
exit(O); 

break;  /*  break  case  F7  */ 


case  F8: 

enter_one_file(); 

strcpy(codes,  "FRQ  2000  HZ,  AMP  10  MV"); 
error  =  change_setrings(codes); 
eiTor_handler  (error,  "lOOUTPUTS"); 
_setvideomode(  _DEFAULTMODE ); 
exit(O); 

break;  /*  break  case  F8  */ 

default: 

control_panel(rpt,  bur,  frq,  amp); 

}  /*  end  switch  */ 

}  /*  end  while  (key  !=  ESC )  */ 

}  I*  end  control_panel()  */ 


/* - - - 

I  Function:collect(numchan,  numbuf,  numtrials,  bsize,  clk_rate) 

.... - - - -  - - ♦/ 

void  collect(int  numchan,  int  numbuf,  int  numtrials,  int  bsize,  int  clk_rate) 
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int  i,  j,  buf_num; 

for  (j  =  0;  j  <  bsize;  j  ++) 

Icxialbuflj]  =  0.; 

for  (i  =  0;  i  <  numtrials;  i++) 

{ 

_settextposition(  29, 1 ); 
sprintf(str_buffer, "Trial  number:  %d.",  i); 

_outtext(  str_buffer ); 

while  ((buf_num  =  sample(numchan,  numbuf,  bsize,  clk_rate))  ==  -1); 
for  (j  =  0;  j  <  bsize;  j  ++) 

localbuf[j]+=((float)buffer[buf_num][j]-2048.)/((float)numtrials*2048.); 

} 

)  7=^  end  collect 


I  Function:  sample(int  numchan,  int  numbuf,  int  bsize,  int  clk_rate) 

. *! 

int  sample(int  numchan,  int  numbuf,  int  bsize,  int  clk_,rate) 

{ 

int  istat,  unit  =  0,  handle,  bhandle,  flag; 
int  i,  j,  k,  hdl[NUMBUF],  channel,  gain; 
int  buf_num; 

char  *driver_name  =  "DT2812$0"; 

BOARDS  board; 

SET_ACQ_PARAMS  sap; 


flag  =  0; 


!*****  initialize  and  reset  the  DT2812  board. 
istat  =  dt_initialize  (driver_name,  &handle); 
if  (istat  !=  0) 

{ 

flag  =  -l; 

_settextposition(  28, 1 ); 

sprintf(str_buffer,”Error  x%x  in  dt_initialize.",  istat); 
_outtext(  str_buffer ); 

) 


/***  create  a  couple  of  buffers  ***/ 
for  (i  =  0;  i  <  numbuf;  i  ++) 

istat  =  dt_create_buffer(unit,AD_SECTION,&bsize,&buffer[i][0]  ,&hdl[i]); 
if  (istat  !=0) 

{ 

flag  =  -l; 

_settextposition(  28, 1 ); 

sprintf(str_buffer,"EiTor  x%x  in  dt_create_buffer.  ",  istat); 

_outtext(  str_buffer ); 
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/*  getchO;  */ 

} 

/*****  create  a  channel/gain  list  *****/ 
channel  =  (numchan  -1)ICHANNEL_SCAN; 
gain  =  0; 

dt_create_cgl  (1,  &channel,  &gain ); 
if  (istat  !=  0) 

{ 

flag  =  -l; 

_settextposition(  28, 1 ); 

sprintf(str_buffer, "Error  x%x  in  dt_create_cgl.  ",  istat); 

_outtext(  str_buffer ); 

) 


/*****  set  board  functions.  *****/ 
board.dma_channel_l  =  0; 
board.interrupt_level  =  0; 
board.board_timeout  =15; 
board.ad_setup_bits  =0x5100; 
board.da_setup_bits  =  0; 

dt_set_board  (unit,  handle,  &board); 
if  (istat  !=  0) 

{ 

flag  =  -l; 

_settextposition(  28, 1 ); 

sprintf(str_buffer,"Error  x%x  in  dt_set_board.  ",  istat); 

_outtext(  str_buffer ); 

} 

/*  setup  the  acquisition  parameters  for  a  burst  adc.  */ 
sap.section  =  AD_SECTION; 

sap.transfer_type  =  0x49;  /*  continuous,  triggered  single  channel  DMA  */ 
sap.clock_rate  =  CLK_RATE; 
dt_set_acq  (unit,  handle,  AD_SECTION,  &sap); 
if  (istat  !=  0) 

{ 

flag  =  -l; 

_settextposition(  28, 1 ); 

sprintf(str_buffer,"Error  x%x  in  dt_set_acq.  ",  istat); 

_outtext(  str_buffer ); 

} 

/*****  start  collection  operation.  *****/ 

istat  =  dt_start_acq  ( unit,  handle,  AD_SECnON); 

if  (istat  !=  0) 

{ 

flag  =  -l; 

_settextposition(  28, 1 ); 

sprintf(str_buffer,"Error  x%x  in  dt_start_acq.  ",  istat); 

_outtext(  str_buffer ); 

} 
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istat  =  dt_wait_buffer(unit,  AD_SECTION,  &bhandle); 
if  (istat  !=  0) 

{ 

flag  =  -l; 

_settextposition(  28, 1 ); 

sprintf(str_buffer, "Error  x%x  in  dt_wait_buffer.  ",  istat); 

_outtext(  str_buffer ); 

} 

buf_num  =  0; 

while  (  bhandle  !=  hdl[buf_num])  buf_num  ++; 

dt_reset_buffer(bhandle); 
if  (istat  !=  0) 

{ 

fiag  =  -l; 

_settextposition(  28, 1 ); 

sprintf(str_buffer,"Error  x%x  in  dt_reset_buffer.  ",  istat); 

_outtext(  str_buffer ); 

/*  getchO;  */ 

) 


!*****  now  stop  the  a/d  conversions  *****/ 
dt_stop_acq  ( unit,  handle,  AD_SECTION ); 
for  (i  =  0;  i  <  numbuf;  i  ++) 

{ 

dt_delete_buffer(unit,  AD.SECTION,  hdl[i]); 
if  (istat  !=  0) 

{ 

flag  =  -l; 

_settextposition(  28, 1 ); 

sprintf(str_buffer,”Error  x%x  in  dt_delete_buffer,  buffer  number  %d.", 

istat,  i); 

_outtext(  str_buffer ); 

} 

} 

!*****  all  done,  terminate  toolkit/driver  communication.  *****/ 

istat  =  determinate  (handle); 

if(istat!=0) 

{ 

flag  =  -l; 

_settextposition(  28,  1 ); 

sprintf(str_buffer, "Error  x%x  in  determinate.  ",  istat); 

_outtext(  str_buffer ); 

) 

if  (flag  =  0)  return  (buf_num);  else  return  (-1); 

)  /*  end  sampleO  */ 


/* - 

I  Function;  enter_data_file0 
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_ *! 

enter_data_fileO 

{ 

int 

error, 

char 

♦codes, 

option, 

destination[50], 
thickness_buffer[3] , 
weight_buffer[  1 0] , 
moisture_buffer[3] , 
bucket_buffer[3], 
compaction_buffer[3] , 
sequence_buffer[3], 
time_buffer[10], 
date_buffer[10]; 

_strtime(  time_buffer ); 

_strdate(  date_buffer ); 

_settextcolor(10); 
for  (j  =  0;  j  <  12;  j++) 

{ 

_settextposition(  19+j,  1 ); 

sprintf(str_buffer,"  "); 

_outtext(  str_buffer ); 

} 

_settextposition(  20,  1  ); 

sprintf(str_buffer, "Enter  soil  type.  (3  chars  only) : "); 

_outtext(  str_buffer ); 
scanf("%s",  file_name); 

_settextposition(  21, 1 ); 

sprintf(str_buffer,"Enter  soil  thickness  in  cm.  (2  chars  only) : "); 

_outtext(  str_buffer ); 
scanf("%s",  thickness_buffer ); 

_settextposition(  22, 1 ); 

sprintf(str_buffer,"Enter  soil  weight  in  kg.  : "); 

_outtext(  str_buffer ); 
scanf("%s",  weight_buffer ); 


_settextposition(  23, 1 ); 

sprintf(str_buffer, "Enter  Moisture.  (1  chars  only) ; "); 

_outtext(  str_buffer ); 
scanf("%s",  moisture_buffer ); 

_settextposition(  24, 1 ); 

sprintf(str_buffer,"Enter  Compaction.  (1  chars  only) : "); 
_outtext(  str_buffer ); 
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scanf("%s",  compaction_buffer  ); 

_settextposition(  25, 1 ); 

sprintf(str_buffer,"Enter  Sequence.  (1  chars  only) : "); 

_outtext(  str_buffer ); 
scanf("%s",  sequence_buffer ); 

_settextposition(  26, 1 ); 

sprintf(str_buffer,"Enter  Bucket  Number.  :  "); 

_outtext(  str_buffer ); 
scanf("%s",  bucket_buffer ); 

strcat(file_name,  thickness_buffer); 
strcat(file_name,  moisture_buffer); 
strcat(file_name,  compaction_buffer); 
strcat(file_nairie,  sequence_buffer); 
strcpy(destination,  "c:\\soiI\\data'\"); 
strcat(destination,  file_name); 
strcat(destination,  ".cmt"); 

_settextposition(  27,  1 ); 

sprintf(str_buffer,"  File  Name  =  %s  ",  destination); 

_outtext(  str_buffer ); 

while(option  !=  'y') 

{ 

_settextcolor(l  1); 

_settextposition(  28, 40 ); 

sprintf(str_buffer,"  Is  the  file  name  correct  (y/n)  "); 

_outtext(  str_buffer ); 
option  =  getcheQ; 
switch(option) 

{ 

case  'n': 

enter_data_file(); 
default : 

break; 

}  /*  end  switch  */ 

}  /*  end  while  */ 


if(  !access(  destination,  0 ) ) 

{ 

_settextcolor(10); 
for(j  =  0;j  <  12;j++) 

{ 

_settextposition(  19+j,  1 ); 
sprintf(str_buffer,"  "); 

_outtext(  str_buffer  ); 

_settextposition(  19+j,  35 ); 
sprintf(str_buffer,"  "); 

_outtext(  str_buffer ); 

) 

_settextposition(  21, 1 ); 
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sprintf(str_buffer,"  %s  ",  destination); 

_outtext(  str_buffer ); 

_settextposition(  22, 1 ); 

sprintf(str_buffer,"File  Alredy  Exists,  Hit  any  key  to  exit. "); 

_outtext(  str_buffer ); 

getchO; 

strcpy(codes,  "FRQ  2000  HZ,  AMP  10  MV"); 
error  =  change_settings(codes); 
error_handler  (error,  "lOOUTPUTS"); 

_setvideomode(  _DEFAULTMODE ); 
exit(0); 

) 


{ 

out  =  fopen(destination,"w"); 

fprintfCout,'^  %s 
fprintf(out,"\n  %s 
fprintf(out,"\n  %s 
fprintf(out,'Vi  %s 
fprintf(out,"\n  %s 
fprintf(out,"\n  %s 
fprintf(out,"\n  %s 
fprintf^out,’^  %s 
fprintf(out,"\n  %s 
fclose(out); 

} 


",file_name); 

",thickness_buffer); 

",weight_buffer); 

",moisture_buffer); 

",compaction_buffer); 

",sequence_buffer); 

",bucket_buffer); 

",nme_buffer); 

",date_buffer); 


_settextcolor(15); 
for  (j  =  0;  j  <  10;j++) 

( 

_settextposition(  19+j,  1  ); 

sprintf(str_buffer,"  "); 

_outtext(  str_buffer ); 

} 

return  (0); 

}  /*  end  enter_data_file0  */ 


/* -  - 

I 

1  Function:  enter_one_file() 

1 

- */ 

enter_one_file0 

{ 

int 

i, 

error, 

char 

♦codes. 
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option, 

destination  [50]; 

_settextcolor(10); 
for(j  =  0;  j  <  12;  j++) 

{ 

_settextposition(  19+j,  1 ); 

sprintf(str_buffer,"  "); 

_outtext(  str_buffer ); 

} 

_settextposition(  20, 1 ); 

sprintf(str_buffer,"Enter  file  name  ( name  &  extennsion ) : "); 

_outtext(  str_buffer ); 
scanf("%s",  file_name); 

strcpy(destination,  "c:\\soiKMata\\"); 
strcat(destination,  file_name); 

_settextposition(  27, 1  ); 

sprintf(str_buffer,"  File  Name  =  %s  ",  destination); 

_outtext(  str_buffer ); 

if(  !access(  destination,  0 ) ) 

{ 

_settextcolor(10); 
for  (j  =  0;  j  <  12;  j++) 

{ 

_settextposition(  19+j,  1 ); 
sprintf(str_buffer,"  "); 

_outtext(  str_buffer ); 

_settextposition(  19+j,  35 ); 
sprintf(str_buffer,"  "); 

_outtext(  str_buffer ); 

) 

_settextposition(  21,1); 
sprintf(str_buffer,"  %s  ",  destination); 

_outtext(  str_buffer ); 

_settextposition(  22, 1 ); 

sprintf(str_buffer,"File  Alredy  Exists.  Hit  any  key  to  exit.  "); 

_outtext(  str_buffer ); 

getchO; 

strcpyCcodes,  "FRQ  2000  HZ,  AMP  10  MV"); 
error  =  change_settings(codes); 
error.handler  (eiror,  "lOOUTPUTS"); 

_setvideomode(  _DEFAULTMODE ); 
exit(O); 

} 

else 

{ 

fp  =  fopen(destination,  "w"); 
colIect(NUMCHAN,  NUMBUF,  NUMTRIALS, 
NUMCHAN+SAMPLES,  CLK_RATE); 

for  (i  =  real_trig;  i<S AMPLES;  i++) 
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[ 

/*  fprintf(fp,  "%d  ",i);  *! 

for  (j  =  0;  j  <  NUMCHAN; }++) 

{ 

fprintf(fp,  "%f  ", 

MAXAMP*localbuf[i*NUMCHAN  +j]); 

} 

fprintf(fp,  ’V); 

} 

fclose(fp); 

strcpyCcodes,  "FRQ  2000  HZ,  AMP  10  MV"); 
error  =  change_settings(codes); 
error.handler  (error.  "lOOUTPUTS"); 

_setvideomode(  _DEFAULTMODE ); 
exit(O); 

} 

)  /*  end  enter_one_file()  */ 


/* . 

I  Function:  error_handler  (error,  routine) 

. */ 

error_handler  (error,  routine) 
int  error, 

char  *routine; 

{ 

char  *estring; 

char  ch; 

if  (error  !=NOERR) 

{ 

printf  ("Error  in  call  to  %s  \n",  routine); 
printf  ("  Error  =  %d  :  %s  \n",  error,  errstr(error)); 
printf  ("Press  <RETURN>  to  continue:  "); 
scanf  ("%c",  &ch); 

) 

}  /*  end  error_handler  */ 


/♦ . . . 

I  Function:  timer() 

I  Purpose :  pause  for  x-seconds  between  taking  data 
1  so  the  signal  can  settle. 

- */ 

timerO 

{ 


ftime(&xtime); 

time  1  =Gong)xtime.  time; 

while(l) 
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{ 

ftime(&xtime); 
time2=(long)xtiine.time; 
if  ( time2-timel  >  wait_time-l) 
break; 

} 

}  /*  end  timer*/ 


Appendices:  Page  A42 


Final  Report 


Appendix  AlO.  Procedure  to  transfer  raw  data  from  soil  acquisition  sytem  to  worksatation. 

After  taking  soil  data  stored  in  c:\soiMata\ . 

From  bmode.brl.uiuc.edu  (128.174.211.58) 


TYPE  at  the  PC 

c:\cd  c:\ncsa2\ftp  -i  128.174.211.3 

username  :  your  name 

password  :  xxxxxx 

ftp>  led  c:\ 

f^>  led  soiI\data 

ftp>  cd  /bugs/soil/raw_data 

ftp>  mput  *.* 

ftp>  quit 


ACTION 

ftp  into  ecstasy  in  interactive  mode 
login  t 

change  into  c:\soiMata\on  PC 

put  data  into  soil  raw  data  directory 

put  all  files  with  the  test...  prefix 

files  from  PC  should  be  over  at  the  SPARC 


The  files  should  now  be  on  ECSTASY  in  the  directory  bugs/soil/raw_data. 
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Appendix  All.  Matlab  program  procsoil.m. 

%  procsoil.m 
% 

%  Produce  time  delay  and  correlation  amplitude  values 
%  from  data  collected  for  US  ARMY  CERL  soil  project. 

% 

%  Version  2.0  January  1996  by  RNC  &  NBS 
% 

%  Saves  data  into  one  big  file  for  processing. 


clear 

soilfile  =  inputCEnter  file  name,  enclosed  in  single  quotes: '); 

thickness  =  eval(soilfile(4:5)); 

moisture  =  eval(soilfile(6)); 

compaction  =  eval(soiIfile(7)); 

sequence  =  eval(soilfile(8)); 

%  check  if  info  is  true! 
dispCthickness  =  ');  disp(thickness); 

dispC'Do  you  want  to  change  this  setting?'); 

ch  =  input(Type  new  thickness  or  -1  to  retain  this  value.'); 

if  (ch  ~=  -1)  thickness  =  ch;  end; 

limitl  =  round(thickness); 

%limitl  =round(1.25*thickness); 
limit2  =  round(5*thickness); . 

open_file  =  [soilfile,  '.prc']; 
open_file(4:5)  =  'xx'; 
fid  =  fopen(  open_file,  'a'); 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

[pk_max_ind,  pk_max_val,  pk_hi_ind,  pk_hi_val,  pk_lo_ind,  pk_lo_val,  pk_med_ind, 
pk_med_val]  =  ... 

peak([soilfile,  '.200',soilfile,  '.201'  soilfile,  '.202',  soilfi 
le,  '.203'],  limitl,  limit2); 

fprintf(fid,  '1000  %l.lf  %1.0f  %1.0f  %1.0f  %1.15f  %1.15f  %1.15f  %1.15f  %1.15f  %1 
.15f  %1.15f  %1.15f«n',  thickness,  moisture,  compaction,  sequence,... 
pk_max_ind/25000,  pk_max_val,  pk_hi_ind/250()0,  pk_hi_val,  pk_med_ind/25000,  pk_m 
ed_val,  pk_lo_ind/25000,  pk_lo_val); 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%%%%%%%%%%%%7o%%%%%%%%%%%%%%%%%%%%%%%%%% 

[pk_max_ind,  pk_max_val,  pk_hi_ind,  pk_hi_val,  pk_lo_ind,  pk_lo_val,  pk_med_ind, 
pk_med_val]  =  ... 

peak([soilfile,  '.204’,soilfile,  ’.205'  soilfile,  '.206',  soilfi 
le,  ’.207'],  limitl,  limit2); 
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fprintf(fid,  '2000  %l.lf  %1.0f  %1.0f  %1.0f  %1.15f  %1.15f  %1.15f  %1.15f  %L15f  %1 
.15f  %1.15f  %1.15f\n',  thickness,  moisture,  compaction,  sequence,... 
pk_max_ind/25000,  pk_max_val,  pk_hi_i  nd/25000,  pk_hi_val,  pk_med_ind/25000,  pk_m 
ed_val,  pk_lo_ind/25000,  pk_lo_val); 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%%%%%%,%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%%%%%%%%%%%%%%%%%%%%%%%7o%%%%%%%%%%%%%%% 

[pk_max_ind,  pk_max_val,  pk_hi_ind,  pk_hi_val,  pk_lo_ind,  pk_lo_val,  pk_med_ind, 
pk_med_val]  = ... 

peak([soilfile,  '.208',soilfile,  '.209'  soilfile,  '.210',  soilfi 
le,  '.211'],  limitl,  limit2); 

fprintf(fid,  '3000  %l.lf  %1.0f  %1.0f  %1.0f  %1.15f  %1.15f  %1.15f  %1.15f  %1.15f  %1 
.15f  %1.15f  %1.15f\n',  thickness,  moisture,  compaction,  sequence,... 
pk_max_ind/25000,  pk_max_val,  pk_hi_ind/25000,  pk_hi_val,  pk_med_ind/25000,  pk_m 
ed_val,  pk_lo_ind/25000,  pk_lo_val); 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

[pk_max_ind,  pk_max_val,  pk_hi_ind,  pk_hi_val,  pk_lo_ind,  pk_lo_val,  pk_med_ind, 
pk_med_val]  =  ... 

peak([soilfile,  '.212', soilfile,  '.213'  soilfile,  '.214',  soilfi 
le,  '.215'],  limitl,  limit2); 

fprintf(fid,  '4000  %  1 . 1  f  %  1  .Of  %  1  .Of  %  1  .Of  %  1 . 1 5f  %  1 . 1 5f  %  1 . 1 5f  %  1 . 1 5f  %  1 . 1 5f  %  1 
.ISf  %1.15f  %1.15f\n',  thickness,  moisture,  compaction,  sequence,... 
pk_max_ind/25000,  pk_max_val,  pk_hi_ind/25000,  pk_hi_val,  pk_med_ind/25000,  pk_m 
ed_val,  pk_lo_ind/25000,  pk_lo_val); 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

[pk_max_ind,  pk_max_val,  pk_hi_ind,  pk_hi_val,  pk_lo_ind,  pk_lo_val,  pk_med_ind, 
pk_med_val]  =  ... 

peak([soilfile. '.216',soilfile, '.217' soilfile, '.218',  soilfi 
le,  '.219'],  limitl,  limit2); 

fprintfffid,  '5000  %l.lf  %1.0f  %1.0f  %1.0f  %r.l5f  %1.15f  %1.15f  %1.15f  %1.15f  %1 
.15f  %1.15f  %1.15f\n',  thickness,  moisture,  compaction,  sequence,... 
pk_max_ind/25000,  pk_max_val,  pk_hi_ind/25000,  pk_hi_val,  pk_med_ind/25000,  pk_m 
ed_val,  pk_lo_ind/25000,  pk_lo_val); 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

[pk_max_ind,  pk_max_val,  pk_hi_ind,  pk_hi_val,  pkjojnd,  pk_lo_val,  pk_med_ind, 
pk_med_val]  =  ... 

peak([soilfile,  '.220',soilfile,  '.221'  soilfile,  '.222',  soilfi 
le,  '.223'],  limitl,  limit2); 

fprintf(fid,  '6000  %l.lf  %1.0f  %1.0f  %1.0f  %1.15f  %1.15f  %1.15f  %1.15f  %1.15f  %1 
.15f  %1.15f  %1.15f\n',  thickness,  moisture,  compaction,  sequence,... 
pk_max_ind/25000,  pk_max_val,  pk_hi_ind/25000,  pk_hi_val,  pk_med_ind/25000,  pk_m 
ed_val,  pk_lo_ind/25000,  pk_lo_val); 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%%%%%%%%7o%%%%%%%%%%%%%%%%%%%%%%%%% 
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%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

[pk_max_ind,  pk_max_val,  pk_hi_ind,  pk_hi_val,  pk_lo_ind,  pk_lo_val,  pk_med_ind, 
pk_med_val]  =  ... 

peak([soilfile,  '.224',soilfile,  '.225'  soilfile,  '.226',  soilfi 
le,  '.227'],  limitl,  limit2); 

fprintf(fid,  '7000  %l.lf  %1.0f  %1.0f  %1.0f  %1.15f  %1.15f  %1.15f  %1.15f  %1 

.15f  %1.15f  %1.15f\n',  thickness,  moisture,  compaction,  sequence,... 
pk_max_ind/25000,  pk_max_val,  pk_hi_ind/25000,  pk_hi_vi,  pk_med_ind/25000,  pk_m 
ed_val,  pk_lo_ind/25000,  pk_lo_val); 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

[pk_max_ind,  pk_max_val,  pk_hi_ind,  pk_hi_val,  pk_lo_ind,  pk_lo_val,  pk_med_ind, 
pk_med_val]  =  ... 

peak([soiirile,  '.228',soiirile,  '.229'  soilfile,  '.230',  soilfi 
le,  '.231'],  limitl,  limit2); 

fprintf(fid,  '8000  %l.lf  %1.0f  %1.0f  %1.0f  %1.15f  %1.15f  %1.15f  %1.15f  %1.15f  %1 
.15f  %1.15f  %1.15f\n',  thickness,  moisture,  compaction,  sequence,... 
pk_max_ind/25000,  pk_max_val,  pk_hi_ind/25000,  pk_hi_val,  pk_med_ind/25000,  pk_m 
ed_val,  pk_lo_ind/25000,  pk_lo_val); 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

[pk_max_ind,  pk_max_val,  pk_hi_ind,  pk_hi_val,  pk_lo_ind,  pk_lo_val,  pk_med_ind, 
pk_med_val]  = ... 

peak([soilfile,  '.232',soilfile,  '.233'  soilfile,  '.234',  soilfi 
le,  '.235'],  limitl,  limit2); 

fprintf(fid,  '9000  %l.lf  %1.0f  %1.0f  %1.0f  %1.15f  %1.15f  %1.15f  %1.15f  %1.15f  %1 
.15f  %1.15f  %1.15f\n',  thickness,  moisture,  compaction,  sequence,... 
pk_max_ind/25000,  pk_max_val,  pk_hi_ind/25000,  pk_hi_val,  pk_med_ind/25000,  pk_m 
ed_val,  pk_lo_ind/25000,  pk_lo_val); 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

fclose(fid); 
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Appendix  A12.  Matlab  program  loadsoil.m. 

function  data  =  loadsoil(filename) 

%  function  data  =  loadsoil(filename) 

% 

%  Load  a  soil  data  file,  regardless  of  extra  characters 

%  at  the  end  of  each  line. 

% 

%  Kate  Frazier  and  Kay  Raum  figured  out  how  to  do  this. 

% 

eval(['fp  =  fopenC".  filename,  "',"r","g");']); 
data  =  fscanf(fp,  '%f,  [4,  inf]); 
data  =  data'; 
fclose(fp); 
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Appendix  A13.  Matlab  program  peak.m. 

function  [pk_max_ind,  pk_max_val,  pk_hi_ind,  pk_hLval,  pk_lo_ind,  pk_lo_val,  pk 
_med_in<i  pk_med_val]  =  peak(files,  limitl,  limit2); 

%  function  [pk_max_ind,  pk_max_val,  pk_hi_ind,  pk_hi_val,  pk_lo_ind,  pk_lo_val, 
pk_med_ind,  pk_med_val]  =  peak(files,  limitl,  limitl); 

% 


file  =  loadsoil(files(l:12)); 
tempi  =file(l:200,l); 
tempi  =  file(l:200,2); 
tempi  =  tempi  -  mean  (tempi); 
tempi  =  tempi  -  mean(temp2); 
cor  =  xcorr(templ,  tempi)'; 
maxes  1  =  maxes(cor); 

file  =  loadsoil(files(13:24)); 
tempi  =  file(l:200,l); 
tempi  =  file( 1:200,2); 
tempi  =  tempi  -  mean  (tempi); 
tempi  =  tempi  -  mean(temp2); 
cor  =  xcorr(templ,  tempi)'; 
maxes!  =  maxes(cor); 

file  =  loadsoil(files(25:36)); 
tempi  =  file(l:200,l); 
tempi  =  file( 1:200,2); 
tempi  =  tempi  -  mean(templ); 
tempi  =  tempi  -  mean  (tempi); 
cor  =  xcorr(templ,  tempi)'; 
maxes3  =  maxes(cor); 

file  =  loadsoil(files(37;48)); 
tempi  =  file(l:200,l); 
tempi  =  file(l:200,2); 
tempi  =  tempi  -  mean(templ); 
tempi  =  tempi  -  mean(temp2); 
cor  =  xcorr(templ,  tempi)'; 
maxes4  =  maxes(cor); 

if  (maxesl(l)  =  0)  maxesl  =  maxesl(2:length(maxesl));  end; 
if  (maxes2(l)  =  0)  maxesl  =  maxes2(2:length(maxes2));  end; 
if  (maxes3(l)  ==  0)  maxes3  =  maxes3(2:length(maxes3));  end; 
if  (maxes4(l)  ==  0)  maxes4  =  maxes4(2:length(maxes4));  end; 

ind  =  zeros(l,401); 
ind(maxesl)  =  ones(size(maxesl)); 
ind(maxes2)  =  ind(maxes2)  +  ones(size(maxes2)); 
ind(maxes3)  =  ind(maxes3)  +  ones(size(maxes3)); 
ind(maxes4)  =  ind(maxes4)  +  ones(size(maxes4)); 

ind  =  ind(200+limitl:200+limit2); 

[maxind,  argmaxind]  =  max(ind); 
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allmaxes  =  fmd(ind  ==  maxind); 


file  =  loadsoil(files(l:12)); 

if  (length(allmaxes)  =  1) 

pk_max_ind  =  limitl  -1  +allmaxes; 
pk_max_val  =  coefl(file,  pk_max_ind); 
else 

pk_max_val  =  -100; 
for  i  =  l:length(allmaxes) 

temp  =  coefl(file,  limitl  -1  +allmaxes(i)); 
if  (temp  >  pk_max_val), 
pk_max_val  =  temp; 
pk_max_ind  =  limitl  -1  +allmaxes(i); 

end; 

end 

end; 

%find  locations  of  peaks 
pk_range  =  limitl  -1  -t-find(ind  >  0); 
values  =  zeros(si2e(pk_range)); 

devaluate  corr  function  at  those  points 
for  i  =  l:length(pk_range), 

values(i)  =  coefl(file,  pk_range(i)); 
end; 

%keep  only  those  points  corresponding  to  significant 
%correlation  values. 

pk_range  =  pk_range(find(values  >  0.  l*max (values))); 

%find  upper  and  lower  bounds  on  tof. 
pk_lo_ind  =  pk_range(l); 
pk_hi_ind  =  pk_range(length(pk_range)); 
pk_med_ind  =  median(pk_range); 

pk_lo_val  =  coefl(file,  pk_lo_ind); 
pk_hi_val  =  coefl(file,  pk_hi_ind); 
pk_med_val  =  coen(file,  pk_med_ind); 


%this  is  the  old  way  of  doing  this. 

%if  (max(ind)  >  1), 

%%  pk_range  =  limitl  -1  -i-find(ind  >  1); 

%  pk_lo_ind  =  pk_ran  ge(  1 ) ; 

%  pk_hi_ind  =  pk_range(length(pk_range)); 

%  pk_med_ind  =  median(pk_range); 

% 

%  pk_lo_val  =  coefl(file,  pk_lo_ind); 

%  pk_hi_val  =  coefl(file,  pk_hi_ind); 

%  pk_med_val  =  coeh  (file,  pk_med_ind); 
%else 

%%  pk_range  =  limitl  -f-rind(ind  >  0); 
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%  pkjojnd  =  pk_range(l); 

%  pk_hi_ind  =  pk_range(length(pk_range)); 

%  pk_med_ind  =  median(pk_range); 

% 

%  pk_lo_val  =  coefl(file,  pkJoJnd); 

%  pk_hi_val  =  coefl(file,  pk_hi_ind); 

%  pk_med_val  =  coeh  (file,  pk_nied_ind); 

%end; 
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Appendix  A14.  Matlab  program  process.m. 


function  process(infile,  outfile) 

eval(['load  infile,  '.prc']); 
data  =  eval(infile); 

temp  =  hist(data(:,l),  (1000:250:10000)); 
frequencies  =  750+250*find(temp  ~=  0); 

fid  =  fopen(outfile,  'w'); 

fprintf(fid,  'freq  speed  1  rhol  attenl  rho_attenl  speed2  rho2  atten2  rho_ 
atten2  speedS  rho3  atten3  rho_atten3  speed4  rho4  atten4  rho_atten4  freql  freq2 
lopel  slope2  slope3  slope4\n'); 
end; 

for  i  =  l:length(frequencies), 
freq  =  frequencies(i); 

%  Identify  thicknesses  present  for  each  frequency. 
thick_ind  =  find(data(:,l)  ==  freq); 
thicknesses  =  data(thick_ind,2); 

%  Most  frequent  TOF,  maximum  corr  if  a  tie. 
tofl  =  data(thick_ind,  6); 

[speedl,  rhol]  =  regression(thicknesses,  tofl); 
speedl  =  0.01/speedl; 
ampl  =  data(thick_ind,  7); 

[attenl,  rho_attenl]  =  regression(thicknesses,  20*logl0(ampl)); 
attenl  = -attenl; 

%  Largest  possible  TOF,  corresponding  corr. 
tofZ  =  data(thick_ind,  8); 

[speed2,  rho2]  =  regression(thicknesses,  tof2); 
speed2  =  0.01/speed2; 
amp2  =  data(thick_ind,  9); 

[atten2,  rho_atten2]  =  regression(thicknesses,  20*logl0(amp2)); 
atten2  =  -atten2; 

%  Median  TOF,  corresponding  corr. 
tof3  =  data(thick_ind,  10); 

[speed3,  rho3]  =  regression(thicknesses,  tof3); 
speed3  =  0.01/speed3; 
amp3  =  data(thick_ind,  11); 

[atten3,  rho_atten3]  =  regression(thicknesses,  20*logl0(amp3)); 
atten3  =  -atten3; 

%  Lowest  possible  TOF,  corresponding  corr. 
tof4  =  data(thick_ind,  12); 

[speed4,  rho4]  =  regression(thicknesses,  tof4); 
speed4  =  0.01/speed4; 
amp4  =  data(thick_ind,  13); 

[atten4,  rho_atten4]  =  regressionfthicknesses,  20*logl0(amp4)); 
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atten4  =  -atten4; 

%  Log  decrement  calculation  for  each  frequency. 
if(i<10), 

filename  =  [infile(l:3),  num2str(max(thicknesses)),  infile(6:8), 

'.20',  num2str(i  -1)]; 
else 

filename  =  [infile(l:3),  num2str(max(thicknesses)),  infile(6:8), 

'.2',  num2str(i  -1)] 
end; 

[fireql,  ffeq2,  slopel,  slope2,  slope3,  slope4]  =  log_dec(filename); 

fprintf(fid,  '%1.0f  %1.5f  %1.5f  %1.5f  %1.5f  %1.5f  %1.5f  %1.5f  %1.5f  %1.5 
f  %1.5f  %1.5f  %1.5f  %1.5f  %1.5f  %1.5f  %1.5f  %1.5f  %1.5f  %1.5f  %1.5f  %1.5f  %1.5f 

V,  ... 

freq,  speed  1,  rhol,  attenl,  rho_attenl,  speed2,  rho2,  atten2,  rh 
o_atten2,  ... 

speed3,  rho3,  atten3,  rho_atten3,  speed4,  rho4,  atten4,  rho_atte 

n4,  ... 

freql,  freq2,  slopel,  slope2,  slope3,  slope4); 

end; 

fclose(fid); 
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Appendix  A15.  Matlab  program  log_dec.m. 

function  [freql,  freq2,  slopel,  slope2,  slopeS,  slope4]  =  log_dec(filename); 

%  Set  system  parameters 
lengthfft  =  512; 
samplerate  =  25000; 

%  Load  data 

orig_frq  =  1000  +  250*eval(filename(ll;12)); 
limitl  =  round((5/orig_frq)*25000)  +2; 
limit2  =  limitl  +  31; 
trace  =  loadsoil(filename); 

%  Keep  only  the  32  samples  immediately  following  the  offset  of  the  exciting  puls 
e. 

%  Subtract  out  offset  value,  and  compute  FFT. 

%  Need  special  processing  to  handle  case  of  upsloping  limit... 
s3  =  trace(limitl:limit2,3); 
s3  =  s3  -  mean(s3); 

S3  =  fft(s3,  lengthfft); 

[x,mx]  =  max(abs(S3(l:lengthfft/2))); 
freql  =  (mx/lengthfft)*samplerate; 

[p,  S]  =  polyfit((l:length(s3)),  log(abs(hilben(s3)))’,  1); 
slopel  = -samplerate*(p(l)/freql); 


s4  =  trace(limitl:limit2,4); 
s4  =  s4  -  mean(s4); 

S4=  fft(s4,  lengthfft); 

[x,mx]  =  max(abs(S4(l:lengthfft/2))); 
freq2  =  (mx/lengthfft)*samplerate; 

[p,  S]  =polyfit((l:length(s4)),  log(abs(hilbert(s4)))',  1); 
slope2  =  -  samplerate*  (p(l)/freq2); 


%  Compute  log  decrement  using  3dB  bandwidths 

band3db_3  =  find(abs(S3(l:lengthfft/2))  >  max(abs(S3(l:lengthfft/2)))/sqrt(2)); 
band3db_4  =  find(abs(S4(l;lengthfftj^))  >  max(abs(S4(l:lengthfft/2))ysqrt(2)); 

BW3  =  ((max(band3db_3)  -  min(band3db_3))/lengthfft)*samplerate; 

BW4  =  ((max(band3db_4)  -  rnin(band3db_4))/lengthfft)*samplerate; 

slope3  =pi*BW3/freql; 
slope4  =  pi*BW4/freq2; 


Appendices:  Page  A53 


Final  Report 


Appendix  A16.  Matlab  program  coef.m. 


%  coef.m 

%  Compute  the  amplitude  of  the  correlation  function  of  the 
%  first  two  columns  of  a  matrix,  evaluated  at  a  given  point. 

function  [value]  =  coef(templ,  delay) 

N  =  200; 

trace  1  =  temp  1(1  :N,1); 
trace!  =  templ(l:N,2); 

trace  1  =  trace  1  -  mean(tracel); 
trace!  =  trace!  -  mean(trace!); 

cor  =  xcorr(tracel,  trace!); 

value  =  cor(N  +delay)/sum(tracel.'^!); 
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Appendix  A17.  Matlab  program  maxes.m. 


function  [ind]  =  maxes(x) 

%  return  indices  of  local  maxima  in  a  data  stream 

z  =  conv(x,  [1  3  1]); 

zl  =  [z  0]; 
z2  =  [0  z]; 
z3  =  zl  -  z2; 

1  =  length(z3); 

z4  =  z3(l:l-l); 
z5  =  z3(2:l); 

z6  =  z4  .*  z5; 

ind  =  find(z6  <=  0); 

ind_temp  =  fmd(z4(ind)  >  z5(ind)); 

ind  =  ind(ind_temp); 

for  i  =  l:length(ind), 
ind(i)  =  ind(i)  - 1 ; 

end; 
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Appendix  A18.  Matlab  program  regression. m. 

function  [slope,  coef]  =  regression(index,  vals) 

[p,  S]  =  poIyfit(index,  vals,  1); 
slope  =  p(l); 

newvals  =  polyval(p,  index,  1); 

coef  =  sum(  (newv^s  -mean(newvals))  .*  (vals  -  mean(vals)))/ ... 
sqrt(sum((newvals  -mean(newvals)).^2)  *  sum((vals  -mean(vals)).'^2)); 
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Appendix  A19.  Mean  acoustic  and  soil  properties. 


Sample 

File 

Number 


Perme¬ 

ability 

10-*cm'' 


adalll 


adall3 


adaWl 


ada211 


ada241 


ada311 


acla341 


adaSll 


cabin 


cabll2 


cabl41 


cab211 


cab241 


cab311 


cab341 


cabSl 1 


drain 


drall4 


dral41 


dral43 


dra211 


dra212 


dra241 


dra311 


dra341 


draSll 


meal  11 


meal  12 


meal41 


mea211 


mea241 


mea3n 


mea341 


meaSll 


lall4 


lallS 


la216 


la217 


la241 


la311 


la341 


la511 


sad  11 


IS 

IS 


0.341 


0.276 


0.200 


0.684 


0.576 


0.393 


0.707 


0.494 


0.172 


0.098 


0.594 


0.364 


0.914 


0.450 


0.897 


0.582 


0.274 


0.163 


0.538 


0.513 


0.499 


0.283 


0.555 


0.505 


0.710 


0.579 


0.142 


0.104 


0.631 


0.391 


0.619 


0.277 


0.903 


0.311 


0.207 


0.230 


0.232 


0.476 


0.495 


0.807 


0.136 


0.753 


0.392 
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Sample 

File 

Number 

Speed 

m/s 

Attn 

Coef 

dB/cm- 

kHz 

Dry 

Soil 

Den 

gm/cc 

Total 

Wet 

Den 

gm/cc 

Soil 

Part 

Den 

gm/cc 

Total 

Solid 

% 

Total 

Pores 

% 

HjO 

Filled 

Pores 

% 

% 

H20 

(vol) 

% 

H20 

(dry) 

% 

H20 

(wet) 

Pore 

Radius 

urn 

Perme¬ 

ability 

10-*cm^ 

sad  12 

164 

0.377 

1.28 

1.30 

2.65 

48.3 

51.7 

2.3 

1.8 

1.8 

6.28 

0.11 

sacl41 

139 

0.363 

1.41 

1.43 

2.65 

53.1 

46.9 

5.5 

2.6 

1.8 

1.8 

6.20 

0.05 

112 

0.443 

1.04 

1.17 

2.65 

39.4 

60.6 

20.5 

12.4 

11.9 

10.6 

6.48 

0.54 

sac212 

JEM 

0.477 

0.92 

1.02 

2.65 

34.8 

65.2 

15.8 

10.3 

11.2 

10.1 

6.79 

1.26 

sac241 

122 

0.687 

1.36 

1.51 

2.65 

51.2 

48.8 

32.3 

15.7 

11.6 

10.4 

6.23 

0.07 

sacSll 

121 

1.24 

1.52 

2.65 

46.9 

53.1 

51.8 

27.5 

22.2 

6.30 

0.14 

sac341 

176 

0.957 

1.62 

1.94 

2.65 

61.3 

38.7 

82.9 

32.0 

19.7 

16.5 

6.10 

0.01 

sacSll 

207 

0.380 

1.35 

1.84 

2.65 

50.9 

49.1 

100.1 

49.1 

36.4 

26.7 

6.23 

0.07 
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